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Abstract 

In this paper we present the computational model underlying the one-way quantum computer 
' which we introduced recently [Phys. Rev. Lett. 86, 5188 (2001)]. The one-way quantum computer 

has the property that any quantum logic network can be simulated on it. Conversely, not all ways of 
quantum information processing that are possible with the one-way quantum computer can be un- 
£N) ' derstood properly in network model terms. We show that the logical depth is, for certain algorithms, 

^ , lower than has so far been known for networks. For example, every quantum circuit in the Clifford 

group can be performed on the one-way quantum computer in a single step. 

O . 

■ 1 Introduction 

Quantum computation has been formulated within different frameworks such as the quantum Turing 
machine |l| or the quantum logic network model 0. In the latter it is particularly easy to establish 
a connection between physics and the processing of quantum information, since the building blocks 
of the quantum logic network -the quantum gates- are unitary transformations generated by suitably 
tailored Hamiltonians. More recently, the capability of projective (von Neumann) measurements to drive 
C3 a quantum computation has been investigated ||-|^). 

In JTo| | we have shown that universal quantum computation can be entirely built on one-qubit mea- 
surements on a certain class of highly entangled multi-qubit states, the cluster states [|ll]. In this scheme, 
a cluster state forms a resource for quantum computation and the set of measurements form the pro- 
gram. This scheme we called the "one-way quantum computer" since the entanglement in a cluster state 
is destroyed by the one-qubit measurements and therefore the cluster state can be used only once. To 
stress the importance of the cluster state for the scheme, we will use in the following the abbreviation 
QC C for "one-way quantum computer" . 

If a quantum logic network is simulated on a QC C , for many algorithms the number of computational 
steps scales more favorably with the input size than the number of steps in the original network does. 
To be specific, circuits which realize transformations in the Clifford group -which is generated by all the 
CNOT-gates, Hadamard-gates and 7r/2-phase shifts- can be performed by a QC C in a single time step, 
i.e. all the measurements to implement such a circuit can be carried out at the same time. Generally, in 
a simulation of a quantum logic network by a one-way quantum computer, the temporal ordering of the 
gates of the network is transformed into a spatial pattern of measurement bases for the individual qubits 
on the resource cluster state. For the temporal ordering of the measurements there is no counterpart in 
the network model. Therefore, the question of complexity of a quantum computation must be possibly 
revisited. 

It should be understood that both the quantum logic network computer and the QC C can simulate 
each other efficiently. The fact that each quantum logic network can be simulated on the QC C has 
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been shown in |T(j]. The converse is also true because a resource cluster state of arbitrary size can be 
created by a quantum logic network of constant logical depth. Furthermore, the subsequent one-qubit 
measurements arc within the set of standard tools employed in the network scheme of computation. In 
this sense, the QC C does not add physical means to quantum computation. However, while the network 
model can describe the means that are used in a computation on a QC C , it cannot describe how they 
have to be used. First, the network for the creation of the cluster state docs not tell anything about the 
computational process itself, since the cluster state is a universal resource. Second, to the description of 
the computational process there belongs the temporal order in which the measurements are performed. 
But the temporal ordering of measurements to simulate quantum gates on the QC C is not pre-imposed 
by the temporal ordering of these gates in the corresponding quantum logic network. For example, in the 
network model two gates cannot be performed in parallel if they do not commute. In the QC c -realization 
they still can if they both belong to the Clifford group. As a consequence, all circuits in the Clifford group 
can be parallelized to logical depth D = 1 on the QC C , as mentioned before. Similarly, all circuits which 
consist of CNOT-gates and either x- or z-rotations with variable angle, can be parallelized to logical 
depth D = 2, independent of the number of logical qubits or gates. The best networks that have been 
found for the described circuits or at least for special cases thereof have logarithmic depth in the number 
of qubits or gates, respectively Jl3| |. 

These observations give rise to two questions: I) "What is the adequate computational model for 
the one-way quantum computer ?" , II) "How can the complexity of computations with the QC C be 
calculated?" . The present paper deals with these two questions. It describes the computational model 
underlying the QC C and provides the tools by which the logical depth of algorithms on the QC C can be 
discussed quantitatively. 

The paper is organized as follows. Section || contains a summary of the QC C , as described in jio]]. In 
Section |3| the terminology required to describe the computational model for the QC C is developed. The 
central part of this paper is Section ^ where the computational model underlying the one-way quantum 
computer is presented. In Section |^, the non-network character of the QC C is illustrated at the hand of 
the temporal complexity for the above mentioned two special classes of circuits. Section |^ is the discussion 
of the results, Section R the conclusion. 



2 Network picture of the QC C 

Before starting the description of the objects relevant for the computational model we would like to give 
a short summary of the one-way quantum computer. In [|lo|| it was proved that the QC C is universal by 
showing that it can simulate any quantum logic network. 

For the one-way quantum computer, the entire resource for the quantum computation is provided 
initially in the form of a specific entangled state -the cluster state jn]]- of a large number of qubits. 
Information is then written onto the cluster, processed, and read out from the cluster by one-particle 
measurements only. The entangled state of the cluster thereby serves as a universal "substrate" for 
any quantum computation. Cluster states can be created efficiently in any system with a quantum 
Ising-typc interaction (at very low temperatures) between two-state particles in a lattice configuration. 
More specifically, to create a cluster state |$)c the qubits on a cluster C are at first all prepared in an 
individual state |+) = l/-\/2~(|0) + |l)) and then brought into a cluster state by switching on the Ising-type 
interaction i/j n t for an appropriately chosen finite time span T. The time evolution operator generated 
by the Ising-type Hamiltonian which takes the initial product state to the cluster state is denoted by S. 
The unitary transformation S on a two-dimensional array of qubits, as used in this paper, has the form 

S = (T^X^ec l°)^°l®l 1 )r-(fc)< 1 l + |0> fe <0|(»|l>„ (fc) <l| 

-in/4 J2 a^a^+a^a^ (1) 



J-J e -iTr/2 o-< fc » e iir/4 ai rW) giTr/4 <T<" (fc)) 

feec 



kec 



where r(k) denotes the site of the right neighbor of qubit k, i.e. the site next following k in the ^-direction 
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and u(k) the site of the upper neighbor of k, i.e. the site next following k in the ^-direction. The 
interaction part of S is generated by an Ising Hamiltonian. In (|l|) there appear additional posterior local 
unitary transformations which have no influence on the entanglement properties of the states generated 
under S. 

The quantum state |$)c, the cluster state of a cluster C of neighbouring qubits provides in advance 
all entanglement that is involved in the subsequent quantum computation. It has been shown fll|| that 
the cluster state |$)c is characterized by a set of eigenvalue equations 

4°) (g) ^'>i$) c = (-in<i>) c , (2) 

a' Gngbh(a) 

where ngbh(a) specifies the sites of all qubits that interact with the qubit at site a G C. The eigenvalues 
are specified by the distribution of the qubits on the lattice and encoded in {k„ G {0, 1}, a = 1, .., ||C|| }. 
They can be altered e.g. by applying phase-flips a z before or after the Ising interaction. For the special 
case of C = Z 2 , that is the case of an infinitely extended cluster, Kk = Vfc G C. The equations (||) 
are central for the described computation scheme. It is important to realize that information processing 
is possible even though the result of every individual measurement in any direction of the Bloch sphere 
is completely random. The reason for the randomness of the measurement results is that the reduced 
density operator for each qubit in the cluster state is |l. While the individual measurement results are 
irrelevant for the computation, the strict correlations between measurement results inferred from (^) are 
what makes the processing of quantum information on the QC C possible. 

Let us for clarity emphasize that in the scheme of the QC C we distinguish between physical cluster 
qubits in C which are measured in the process of computation, and the logical qubits. The logical qubits 
constitute the quantum information being processed while the cluster qubits in the initial cluster state 
form an entanglement resource. Measurements of their individual one-qubit state drive the computation. 

To process quantum information with this cluster, it suffices to measure its particles in a certain 
order and in a certain basis. Quantum information is thereby propagated through the cluster and 
processed. Measurements of cr z -observables effectively remove the respective lattice qubits from the 
cluster. Measurements of a x are used for "wires" i.e. to propagate logical quantum bits through the 
cluster, and for the CNOT-gate between two logical qubits. Observables of the form cos(tyj) a x + sm((p) a y 
are measured to realize arbitrary rotations of logical qubits. Here, the angle tp specifies the measurement 
direction. For the one-qubit rotations, the basis in which a certain qubit is measured depends on the 
results of preceding measurements. This introduces a temporal ordering in which the measurements have 
to be performed. The processing is finished once all qubits except a last one on each wire have been 
measured. The remaining unmeasured qubits form the quantum register which is now ready for the 
readout. At this point, the results of previous measurements determine in which basis these "output" 
qubits need to be measured for the final readout, or if the readout measurements are in the a x -, a y - 
or (T^-cigcnbasis, how the readout measurements have to be interpreted. Without loss of generality, we 
assume in this paper that the readout measurements are performed in the cr 2 -eigenbasis. 

For illustration and later reference we review two points of the universality proof for the QC C . First, 
the realization of the arbitrary one-qubit rotation and of the CNOT-gate as the elements of the universal 
set of gates. And second, the effect of the randomness of the individual measurement results and how to 
account for them. 

An arbitrary rotation Un t G SU(2) can be achieved in a chain of 5 qubits. Consider a rotation in its 
Euler representation 

u Ro t{UiX) = u x (()umu x (0, (3) 

where the rotations about the x- and z-axis arc 



U x {a) = 
U z {a) = 
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Initially, the first qubit is in some state |V'in)i which is to be rotated, and the other qubits are in |+). After 
the 5 qubits are entangled by the time evolution operator S generated by the Ising-typc Hamiltonian, 
the state (V'm) can be rotated by measuring qubits 1 to 4. At the same time, the state is also swapped 
to site 5. The qubits 1 ... 4 are measured in appropriately chosen bases, viz. 

Bi( Vi „) = | ^ , ^ 1 (5) 

whereby the measurement outcomes Sj € {0, 1} for j = 1 . . .4 are obtained. Here, Sj = means that 
qubit j is projected into the first state of Bj((pj tJaeaa ). In (0) the basis states of all possible measurement 
bases lie on the equator of the Bloch sphere, i.e. on the intersection of the Bloch sphere with the x-y-plane. 
Therefore, the measurement basis for qubit j can be specified by a single parameter, the measurement 
angle <pj imea s- The measurement direction of qubit j is the vector on the Bloch sphere which corresponds 
to the first state in the measurement basis Bj(<pj, m eas)- Thus, the measurement angle ^meas is equal 
to the angle between the measurement direction at qubit j and the positive x-axis. For all of the so far 
constructed gates, the cluster qubits are either -if they are not required for the realization of the circuit- 
measured in cr z , or -if they are required- measured in some measurement direction in the x-y-plane. 
In summary, the procedure to implement an arbitrary rotation »?, C), specified by its Euler angles 

f,77,C, is this: 

1 . measure qubit 1 in B\ (0) 

2. measure qubit 2 in B 2 (-£ (-l) sl+ <') 

3. measure qubit 3 in B3 l—r) (— 1) S2+K 2^ ^ 

4. measure qubit 4 in B 4 (-£ (_l) s i+ s 3+ K i./+ K ^ 

If the 5-qubit cluster state in Fig. [l]is created from a product state of all qubits in |+) via the interaction 
S of eq. (Q), we have the specific values j = 0, n' 2 = n' a = « 4 = k' 5 = 1. Please note that in (^|) we 
used the set {k^ j, k' 2 , k' 3 , k' 4 , k' 5 q } instead of {ki,1 = 1,..,5} to specify a cluster state 1 0) 5 on a chain 
of 5 qubits. We will use primed k, whenever we describe a cluster state \(j)}c N on some cluster Cn 
which consists only of those cluster qubits that are necessary for the implementation of some circuit. 
Further, the k! associated with the input qubit 1 we denote by k[ j instead of k[, and the k' associated 
with the output qubit 5 we denote by k' 5 instead of k' 5 . The reason for this notation will be discussed 



in Section 4.1. The unprimed k, {/t n }, are reserved for the description of the cluster state \<p)c on the 
whole cluster C. 

Via the procedure (0) the rotation U' Rot is realized: 

u'bAL v, = u^ Rot u Rot (£, n, C). (7) 

Therein, the random byproduct operator has the form 

„ S2 + S4+k' 2 +k' 4 S 1 +S 3 +k' 1 j+k' 3 +k' 5 , . 

U\Z,Rot = Ox o z . (8j 

It can be corrected for at the end of the computation, as explained below. 

A CNOT-gate can be implemented on a cluster state of 15 qubits, as shown in Fig. [j]. All measurements 
can be performed simultaneously. Depending on the measurement results, the following gate is thereby 
realized: 

U'cnot = Us,cnot CNOT( Cl t). (9) 
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Therein the byproduct operator Uy,,cnot has the form 

tt ( C ) 7 - C) W 7 *" ( C ) 7 = C> t*) 7 '" vv, 

Uy,,cnot = Ox o x ' <j z ' o z , with 

7x = «2 + «3 + S5 + «6 + «2 + K 3 + K 5 + K 6 

7* = S2 + S3 + S S + S10 + Sl2 + S14 + K 2 + K 3 + K 8 + K 10 + K 12 + K 14 ( 10 ) 



(c) 

7z = S1 + S3 + S4 + S5 + S8 + S9 + S11+ 

7« = S9 + S11 + S 13 +K9 j + k' u +«' 13 +«' 150 . 



Please note the constant offset for 7^' and further that for both the general rotation and the CNOT-gate 
the byproduct operators depend only on the combinations Sk + n' k if k is not an output qubit of the gate. 



The latter is a general feature of the relevant qubits, as will be shown in Section 3.4.1. 

The randomness of the measurement results does not jeopardize the function of the circuit. Depending 
on the measurement results, extra rotations o x and a z act on the output qubit of every implemented 
gate, as in (j?j), for example. By use of the propagation relation for general one-qubit rotations 

U R (^C)o s z oi = ata^ U R ((-ir^(~iy' V ,(-irO, (11) 
and the one for the CNOT-gate, as in |l2| , 

CNOT(c,t)a^ St ai c ~> Sc a^ S,t a x c ^ = a^ +St a^ <+< CNOT(c,i), (12) 

these extra rotations can be propagated through the network to act upon the output state. By this we 
mean that the extra random rotations need not be corrected for after a gate. Instead, one just needs to 
keep track of them and delay correction until the end of the computation. Then, the extra rotations can 
be accounted for by properly interpreting the er z -readout measurement results. 

The propagation relations ( |TT| ) for the arbitrary rotation and ( |l2"| ) for the CNOT-gate differ with 
respect to which of the two unitary transformations -the gate or the byproduct operator t/s- is modified 
on the right hand side of ( [TT] ) and (12). In the case of the propagation of a byproduct operator through a 



rotation (|TTJ), the gate is changed and the byproduct operator remains unchanged. It passes just through. 
Conversely, in the case of the propagation of a byproduct operator Us through a CNOT-gate ([l2]), it is 
the byproduct operator which is modified and the gate remains unchanged. 

The measurement bases B(p) and B{— p) in (g) coincide for angles ip = and for ip = ±ir/2. For 
ip = the measurement basis B((p) is the eigenbasis of o x , and for p = ±7r/2 the measurement basis 
B(tp) is the eigenbasis of a y . In these cases, the choice of the measurement basis is not influenced by 
the results of measurements at other qubits. Therefore, rotations whose Eulcr angles £, rj, £ are in the 
set {0,±7r/2} can be realized simultaneously in the first round of measurements, that is no other cluster 
qubits need to be measured before. Among these rotations are the Hadamard gate and the ir/2 phase 
shift. As displayed in Fig. |l|, the Hadamard gate and the 7r/2-phase shift are both realized by performing 
a pattern of a x - and cry-measurements on the cluster C. The byproduct operators which are thereby 
created are 

S1+S3 + S4 + KJ j+Kj+Kj S2+S3+K2 + 4+ K 5 O 

Uz,h = o x a z 

u T,,U z {k/2) — O x O z 

Owing to the special Eulcr angles for the Hadamard- and the 7r/2-phase gate, the propagation relations 
for these rotations can also be written in a form resembling the propagation relation ( |l2| ) for the CNOT- 
gate 

// o, - a. - = a x s *a z s *H, . . 

U z {n/2)a x ''a,'' = o x s * a z s *+ s * U z (n /2). 1 > 
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Figure 1: Realization of elementary quantum gates on the QC C . Each square represents a 
lattice qubit. The squares in the extreme left column marked with white circles denote the 
input qubits, they on the right-most column denote the output qubits. Note that blank 
squares can represent either a z measurements or empty lattice sites. 



Under propagation -via the propagation relations (|Tl]), ( |l2| ) and ([14])- the byproduct operators resulting 
from the implementation of the universal gates generate a subgroup of the group 

i/ looal = h,crg\(T®,i= l..n, and products thereofj (15) 

of all possible byproduct operators. £/ local c SU(2)® n is a subset of the set of all multi-local unitary 
operations. Hence, it can be compensated for at the end of the computation by a local change of the 
measurement bases. 

To summarize, any quantum logic network can be simulated on a one-way quantum computer. A set 
of universal gates can be realized by one-qubit measurements and the gates can be combined to circuits. 
Due to the randomness of the results of the individual measurements, unwanted byproduct operators 
are introduced. These byproduct operators can be accounted for by adapting measurement directions 
throughout the process of computation. In this way, a subset of qubits on the cluster C is prepared as 
the output register. The quantum state on this subset of qubits equals that of the quantum register of 
the simulated network up to the action of an accumulated byproduct operator. The byproduct operator 
determines how the measurements on the output register are to be interpreted. 

3 Beyond the network picture 

In the previous section we have described the QC C in a network terminology, which has been useful to 
prove the universality of the scheme. On the other hand, the cluster qubits do not have to be measured 
in the order prescribed by the order of the gates in the corresponding network. This observation indicates 
that the network picture does not describe the QC C in every respect. 
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3.1 The sets Q t of simultaneously measurable qubits 



The cluster qubits which we have chosen to take the role of the readout register, for example, are just 
qubits like any other cluster qubits. It turns out that, in a more efficient way of running the QC C , the 
"readout" qubits are not the last ones to be measured but among the first. It is advantageous to forget 
about the network altogether and to view the QC C as a set of one-qubit measurements on a resource 
quantum state, the cluster state. These measurements have to be performed in a certain order and in a 
certain basis. The classical information of how to measure subsequent qubits must all be contained in 
the results of the already performed measurements. Similarly, the final result of the computation must 
be contained in all the measurement outcomes together. 

In the following we will adopt the strategy that every cluster qubit is measured at the earliest possible 
time. This means that each qubit is measured as soon as the required measurement results from other 
qubits which determine its measurement basis are known. Let us denote by Q t the set of qubits which 
can be measured at the same time in the measurement round t. So, how can the sets Q t be determined? 
Qo is the set of qubits which are measured in the first round. These are all the qubits whose observables 
a x , a v or a z are measured. The measurement bases for these qubits do not depend on the results of any 
previous measurements. To determine the subsequent set Qi, one looks at which qubits can be measured 
with the knowledge of the measurement results from the qubits in Qo- Next, one looks which qubits can 
be measured with the measurement results from the qubits in Qo and Qi known. These qubits form the 
set Qi- In this manner one proceeds until the whole cluster C is divided into disjoint subsets Qt- 

As will become clear later, it is useful to introduce the sets Q^ of yet-to-be measured qubits. More 
precisely, Q^ is the set of qubits which remain to be measured after measurement round No. t — 1, 



Mathematically, the sets Qt are derived from a strict partial ordering in C. The strict partial ordering, 
in turn, is generated by forward cones which are explained in the next section. 

3.2 The forward- and backward cones 

Be g a gate in the network Af to be simulated and k G C(g) a cluster qubit that belongs to the imple- 
mentation of g. Further, be O, A and f2 three vertical cuts through the network Af. A vertical cut is 
such that it intersects each qubit line in a network only once and that it does not intersect gates. O 
intersects M just after the gate g, i.e. the byproduct operator U^,g caused by the implementation of g, 
as given in (g) and (0), is located on O. Note that the byproduct operators generated on O depend on 
the measurement results obtained in course of the gate implementation via (Uk) Sk ■ A intersects Adjust 
be fore the input, i.e. an operator propagated to A acts on the input register of Af, and O intersects Af 
just before the output such that an operator propagated to il acts on the output register of Af. We can 
now define the forward- and backward cones of the cluster qubits k G C. 

Definition 1 The forward cone fc(fc) of a cluster qubit k G C is the set of all those cluster qubits j G Q^ 
whose measurement basis B((pj^ mea , s ) depends on the result Sk of the measurement of qubit k after the 
byproduct operator (Uk)' Sk is propagated from O to fi. 

Definition 2 The backward cone hc{k) of a cluster qubit k G C is the set of all those cluster qubits 
j G Q^ whose measurement basis £>(</?j.mcas) depends on the result s^ of the measurement of qubit k 
after the byproduct operator (Uk) Sk is propagated from O to A. 

It will turn out that only the backward cones of the qubits k G Qo constitute part of the information 
specifying an algorithm on the QC C , but nevertheless all the backward- and forward cones are important 




(16) 
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objects in the scheme. Either of the sets, the set of the forward- and that of the backward cones, separately 
contains the full information of the temporal structure of a computation on the QC C . 

Let us examine the definitions [l] and || for a particular example, the general one-qubit rotation (|J) as 
implemented by the procedure ds|) modulo a byproduct operator Us, Rot as given in (||). The measurement 
result si of qubit 1 (cf. Fig. modifies the measurement angle of qubit 2, which is responsible for 
implementing an x-rotation U x {^), by a factor (— l) Sl . Further, it causes a byproduct operator (cr z ) Sl at 
O. If this byproduct operator is propagated forward from O to fl it has no effect on qubit 2, because 
qubit 2 is behind O. The dependence on si of the basis in which qubit 2 has to be measured persists 
and thus qubit 2 is in the forward cone of qubit 1, 2 € fc(l). The situation is different if the byproduct 
operator (er z ) si is propagated backwards from O to A: via the propagation relation ( fi]]) the Euler angle 
£ is modified by a factor (— 1) 1 which has to be accounted for by multiplying the measurement angle 
V2,meas by a factor (— l) Sl , too. Thus, the factor (— l)* 1 modifies the measurement angle </J2,meas twice, 
once via the procedure (^|) and once in backward propagation, and there is no net effect. Qubit 2 is not 
in the backward cone of qubit 1, 2 ^ bc(l). 

What does it mean that a cluster qubit j is in the forward cone of another cluster qubit k, j G 
fc(fc)? According to the definition, a byproduct operator created via the measurement at cluster qubit k 
influences the measurement angle y^meas at cluster qubit j. To determine the measurement angle at j 
one must thus wait for the measurement result at k to see what the byproduct operator created randomly 
by the measurement at k is. If j G fc(fc), the measurement at qubit j is performed later than that at 
qubit k. This we denote by k -< j 

j G fc(fc) => k -< j. (17) 

Please note that the converse of (JTt]) is not true. If k -< j holds, still j G fc(fc) may not. This can be easily 
verified for the example of a general rotation (||). There, according to the procedure for implementing 
such a rotation described in Section |^, the result of the measurement of qubit 1 enters into in which basis 
qubit 2 has to be measured. Hence, 2 <G fc(l). By (|l7j), 1 -< 2 which means that the measurement at 
qubit 2 has to wait for the result of the measurement on qubit 1. Similarly, the measurement result on 
qubit 2 enters in the choice of the measurement basis for the measurement on qubit 3. 3 € fc(2) and thus 
2 -i 3. Then 1^3 also holds as shown below in ( jig ) , but 3 ^ fc(l), since the measurement result on 
qubit 1 does not influence the choice of the measurement basis for the measurement on qubit 3. 

The relation is a strict partial ordering. Suppose, that besides k -< j, for another cluster qubit 
I one had I G fc(j') and thus j -< I. This would mean that the measurement at I must wait for the 
measurement at j, which itself had to wait for the measurement at k. Thus, the measurement at I also 
had to wait for the measurement at k. Therefore the relation "-<" is transitive, 

k<3 A j -<l — > k<l. (18) 

Further, a measurement to implement a gate cannot and does not need to wait for its own result. 
Therefore the relation is anti-reflexive, 

-dj C : j - ./. (19) 

Let us now cast the procedure to construct the sets of simultaneously measured qubits given above in 
more precise terms. Be Q t C C the set of cluster qubits measured in measurement round t, and C C 
the set of qubits which are to be measured in the measurement round t and all subsequent rounds, as 
defined in ([if]). Then, Qo is the set of qubits which are measured in the first round. These are the qubits 
of which the observablcs a x , a y or a z arc measured, so that the measurement bases arc not influenced by 
other measurement results. Further, = C. Now, the sequence of sets Qt can be constructed using 
the following recursion relation 

Qt = { q eQ {t) h3 P eQ {t} : P ^q} (w) 
Q( t+1 > = Q (t) \Q t . 
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All those qubits which have no precursors in some remaining set and thus do not have to wait for 
results of measurements of qubits in Q^' are taken out of this set to form Q t . The recursion proceeds 
until Q(*max+1) = f or 

some maximal value i ma x of t. 

Can it happen that the recursion does not terminate? That were the case if for a number m of qubits 
jij •••! 3 m € C formed a cycle ji -< j% -< ... -< j m -< ji- Then, none of the qubits ji, .., j n could ever taken 
out of the set. However, by transitivity ([l8]) we then had ji -< j\ which contradicts anti-reflexivity ([To|) . 
Hence, such a situation cannot occur. 

Let us at the end of this section define the forward- and backward cones fc(<7), bc(<?) of the gates g. 
In eq. ( |io| ) we have seen that the byproduct operator caused by the implementation of a CNOT-gate 
contains a constant contribution Uq{CNOT) = ai°\ This contribution to the byproduct operator does 
not depend on any local variables such as the measurement results and is thus attributed to the gate as a 
whole. This byproduct operator is of the same form as those depending on the individual measurement 
results and can influence measurement angles when being propagated forward or backward. Thus we 
define the forward- and backward cones of gates, in analogy to those of the cluster qubits k £ C, as 
follows: 

The forward cone fc(g) of a gate g £ J\f is the set of all those cluster qubits j £ of which the 
measurement basis B(tpj imeas ) is modified if the byproduct operator Uo(g) is propagated forward from O 
to 0. 

The backward cone bc(<7) of a gate g £ Af is the set of all those cluster qubits j £ Q^> of which the 
measurement basis B{ipj^ meas ) is modified if the byproduct operator Uo(g) is propagated backward from 
O to A. 

The forward- and backward cones of gates contain some information about the temporal structure of 
an algorithm on the QC C , but -in contrast to the backward- and forward cones of the cluster qubits- 
not all. Also they do not form part of the information representing a quantum algorithm on the QC C , 
they will be absorbed into the algorithm angles. Their role for the description of a computation on the 
QC C is a technical one. 



3.3 The algorithm- and measurement angles 

There are three different types of angles involved in the described scheme of quantum computation of 
which the most prominent are the algorithm angles and the measurement angles. 

The algorithm angles {vj,algoi 3 & Q^} are part of the information that specifies an algorithm on the 
QC C . They are derived from the network angles {<Pj : q\ n ,j £ Q^}, i.e. the Euler angles of the one-qubit 
rotations in the quantum logic network. Further, the algorithm angles depend on the set \Kk,k £ C} 
characterizing the cluster state \4>)c in @), and on special properties of the measurement pattern. We 
see that the network angles are absorbed into the algorithm angles. They do not constitute part of the 
information specifying a QC c -algorithm. 

As described before, the process of computation with the QC C comprises several measurement rounds. 
The first round, in wich the qubits in the set Qo are measured, is somewhat different from the following 
rounds. Therein, all gates of the circuit that belong to the Clifford group are implemented at the same 
time, no matter where they are located in the corresponding quantum logic network and in which step 
they would be carried out there. This results in byproduct operators scattered all over the place. These 



byproduct operators are, according to the scheme described in Section 4.2, propagated backwards. To 
account for the effect that the byproduct operators have on the algorithm angles, these angles have to 
updated to the modified algorithm angles { algo ' ^ Q^}- The modified algorithm angles <// algo are 
calculated from the respective algorithm angles <Pj A \g and the results obtained in the first measurement 
round j Sfe , k £ Qq } . In the subsequent measurement rounds no further update of the modified algorithm 
angles occurs. Finally, each qubit j £ Q t C is measured in some measurement round t in the basis 
i3((/?j jm eas) where <fj,meas denotes the measurement angle of qubit j. The measurement angle (^meas of 
a qubit j £ Q t is calculated from the modified algorithm angle, ip'. al and the results {sfc, k £ Ui=o Qi] 
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of the so far obtained measurements. 

Before a quantum algorithm is run on the QC C , the algorithm angles are determined from the cluster 
and the properties of the algorithm. During runtime of the QC C , in the first measurement round (t = 0), 
the algorithm angles arc replaced by the modified algorithm angles, i.e. only the latter are kept while the 
former are erased. Then, in the measurement round t a qubit j € Qt is measured in the basis determined 
by the measurement angle tpj^aeaa- After the measurement of qubit j both <//. algo and y^meas can be 
erased. 

Now there arises the question of how the measurement angles of the actual measurements are cal- 



culated from the results of previous measurements. This question will be answered in Section 4.2. The 



question which interests us most, of course, is: "How can the final result of the computation be deter- 
mined from all the measurement outcomes?" It will turn out that the answers to both questions are very 
much related. 



3.4 Quantities for the processing of the measurement results 
3.4.1 The information vector 



Let us again -as in Section 3.1- make the change of the viewpoint from an explanation of the QC C within 
the network model to a more suitable description, but now focusing on the quantities required to describe 
the processing of the measurement results. 

Suppose in the simulation of a quantum logic network J\f on the QC C in a network manner -i.e. 
measuring the "readout" qubits at last- the processing has reached the stage where all but those cluster 
qubits have been measured which form the output register. 

The accumulated byproduct operator Us to act upon the logical output qubits 1, ..,n is known. It 
has the form 

n 

^=n (21) 

i=l 

where Xi, zi G {0, 1} for 1 < i < n. Let us now label the unmeasured qubits on the cluster in the same 
way as the readout qubits on the quantum logic network are labelled. 

The "readout" qubits on the cluster are, at this point, in a state [/si out), where |out) is the output 
state of the corresponding quantum logic network. In the network picture, the computation is completed 
by measuring each qubit in the <7 z -eigenbasis, thereby obtaining the measurement results {sj}, say. In 
the QC C scheme, one measures the state C/s|out) directly, whereby outcomes {si} are obtained and the 

"readout" qubits are projected into the state \Ai) — 11™= i 1+ ^~*2 E/j:|out). Depending on the known 
byproduct operator Us, the set of measurement results {s,} in general has a different interpretation from 
what the network readout {s^} would have. The measurement basis is the same. From ( ^l| ) one obtains 

~2~ 



\m) = niu i+ <-y"*' Efcioot) = u s (ui niu 1+( T g - ^) m> (22) 



From ( |22] ) we see that a <7 z -measurement on the state Us |out) with results {s^ represents the same 
algorithmic output as a CTz-measurement of the state |out) with the results {s^}, where the sets {si} and 
{s-} are related by 

s[ = S{ + Xj mod 2. (23) 

The set {s^} represents the result of the computation. It can be calculated from the results {si} of the 
admeasurements on the "readout" cluster qubits, and the values {xi} which can be extracted from the 
byproduct operator Us- We see that both the contribution of the byproduct operator and the result of the 
measurement on the "readout" qubits of the cluster enter expression (^3|) in the same way. Indeed, there 
is no need to distinguish between these two contributions. On the level of the byproduct operators, the 
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readout measurement result is translated into an additional contribution to the accumulated byproduct 
operator. Both contributions to the such extended byproduct operator 



Uzr = UzUr, with U R = Y[ (c 



(24) 



stem from random measurement results. It is just that the contributions which constitute Us must 
be propagated forward from where they originated and the additional contributions from the readout 
measurements must be propagated one step backwards. Both the forward and the backward propagated 
contributions to U-^r are propagated to the same location in the network. Forward and backward propa- 
gation are closely related. In fact, as the propagation relations ([H]), (|l^) and (jlj) are their own inverse, 
the rules are the same for both directions. The distinguished role of the readout qubits is only a remnant 
of the interpretation of the QC C as a quantum logic network. A more adequate description will have the 
consequence that the cluster qubits on the "output register" , for instance, will be measured during the 
initialization of the QC C such that they are removed from the entangled quantum state even before the 
main part of the computation starts. 

We now define the the information vector I, a 2n-component binary vector which is a function of the 
quantities {xi,Zi} and the results {s^} of the measurements on the cluster output register. Together, 
these quantities determine the extended byproduct operator U^r via (21), and([2^). 



Definition 3 The information vector I is given by 



( \ 



I, 

V J 



, with I x 



( X\ + Si \ 

x 2 + s 2 



\ *En "1" Sn J 



( ^ \ 

Z2 



(25) 



As can be seen from ( |23[ ) and (|25|), I x is a possible result of a re adout measurement in a corresponding 
quantum logic network. I z is redundant. However, in Section 4J the flow quantity I(t), the information 
flow vector, will be defined for which I(t m ax) = I, with i max the index of the final computational step. 
For t < tmax, in 1(0 both the z-part l z (t) and the x-part I x (t) are required to determine the bases for 
the one-qubit measurements in Qt+i- As I z (i) is of equal importance as l x (t) throughout the process of 
computation we keep I z in the definition of I as well. 

The set of possible information vectors I forms a 2n dimensional vector space over V. Let us 
consider the group £/ local of all possible extended byproduct operators U-^r. If wc divide out the normal 
divisor {±1} of W ocal , the resulting factor group U = U lacal / {±1} is isomorphic to V. From the viewpoint 
of physics, dividing out the normal divisor {±1} means that we ignore a global phase. The isomorphism 
X which maps an I € V to the corresponding Usr € U is given by 



1: V3l 



n 



rW 



rW 



(26) 



where [I x ]i and [I z ] i are the respective components of l x and I z . The component-wise addition of vectors 
in V corresponds, via the isomorphism X, to the multiplication of byproduct operators modulo a phase 
factor {±1}. The procedure to implement this product is to first use the operator product, then bring the 
factors into normal order according to (^6|) and finally drop the phase. Multiplication of vectors I G V 
with the scalars 0,1 corresponds to raising the byproduct operators Usr € U to the respective powers. 
One may switch between the two pictures via the isomorphism (|2^). The algebraic structures involved 
will be more apparent in the representation using the information vector 1 = 1 {{xi, Zi, Sj}) than in the 
formulation of the operator U^r. 
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Now that we have defined the information vector I in (g5|) and have seen that the result of the 
computation can be directly read off from the rr-part of I, we would like to find out how I depends on the 
measurement outcomes {sfc} and the set {nk} of binary numbers that determine the cluster state \<f>)c 
in (||). This task is left until Section |4.l| . Before we can accomplish it we need some further definitions. 
It will turn out that the information vector I can be written as a linear combination of the byproduct 
images which are explained next. 



3.4.2 The byproduct images 

Be fi the "cut" through a network Af which intersects the qubit lines just before its output. This is the 
cut at which the extended byproduct operator Usr is accumulated. Consider a qubit k on the cluster C 
which is measured in the course of computation. Depending on the result of the measurement on qubit 
fc, a byproduct operator (Uk)" k is introduced in Af at the location of the logical output qubits of the gate 
for whose implementation the cluster qubit k was measured. This byproduct operator Uk can -by using 
the propagation relations (p"i"|), ( |l2| ) and (|Ti|)- propagated from where it occurred to the cut fl. There 
it appears as the forward propagated byproduct operator Uk\n- Now we can define the byproduct image 
Ffc of a cluster qubit k e C. 

Definition 4 Each cluster qubit k e C has a byproduct image which is the vector that corresponds 
via the isomorphism X" 1 ( |26| ) to the forward propagated byproduct operator Uk\a, 

F k =I- 1 (U k \ n ). (27) 



In the definition (27) of the byproduct image F& it is mentioned only implicitly that the image is evaluated 
on the cut O. Later in the discussion it will become apparent that we could evaluate the byproduct image 
on every vertical cut O. Sometimes, if we compare to other vertical cuts, we will explicitly write Ffc|n 
for P fc . 

The set of byproduct images {Fj., k G C} is an important quantity for the scheme. It represents part 
of the information which is needed to run a quantum algorithm with the QC C . 

In eq. (|Io|) the byproduct operator for the CNOT-gate as realized according to Fig. |l| is given. This 
byproduct operator contains a constant contribution Uq(CNOT) = ai c \ As Uq docs not depend on any 
local variables, neither on {s/c} nor on {k^}, it makes no sense to attribute it to any of the cluster qubits 
that were measured to realize the gate. Instead, it is attributed to the part of the measurement pattern 
that implements the gate as a whole, or -for simplicity- to the gate itself. For any gate g, Uo(g) can 
be propagated forward to the cut fl to act upon the "readout" qubits. There it appears as the forward 
propagated byproduct operator Uo(g)\n- In analogy to the byproduct images of the cluster qubits, we 
can now define the byproduct images of the gates g of the quantum logic network that is simulated on 
the QC C . For any such gate g the byproduct image F g is the vector that corresponds to Uo(g)\n via 

F fl =2^0(0)10). (28) 

Please note that in contrast to the byproduct images F^ of cluster qubits k 6 C the byproduct images F g 
of gates do not form a separate part of the information specifying a quantum algorithm on the QC C . They 



will be absorbed into the initialization value Ij n it of the information flow vector defined in Section 4.1 
and they are thus only a convenient tool in the derivation of the computational model. 

Via I -1 we map the multiplication of byproduct operators, i.e. their accumulation, onto addition 
modulo 2 on the level of the vectors in V. Now there arises the question whether other operations on the 
byproduct operators could be expressed in terms of the corresponding vectors, too. Specifically, one may 
ask how the byproduct operator propagation looks like on the level of the I E V. 



3.4.3 The propagation matrices 

The answer to this question is that on the level of the vector quantities in V propagation is described 
by multiplication with certain 2n x 2n-matrices C . Consider two cuts 0\ and C 2 through a network 
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which intersect each qubit line only once. Further, be the two cuts such that they do not intersect each 
other and that 0\ is earlier than 2 - The part of the quantum logic network between 0\ and 2 is 
denoted by Afo 1 ^o 2 - Be lk\oi and lk\o 2 the vectors describing a byproduct operator resulting from the 
measurement of qubit k, propagated to the cuts 0\ and 2 , respectively. Then we have 



)I 



(29) 



To any quantum logic network A/ a matrix CV can be assigned. For a network A/2 ° A/i composed of two 
subnetworks A/i and A/2 (of which A/i is carried out first) the propagation matrix is equal to the product 
of the propagation matrices of the subnetworks 



C{N 2 oN l ) = C{M 2 )C{N 1 ). 



(30) 



Because of property ( p0|) we only need to find the propagation matrices for the general one-qubit rotations, 
the CNOT-, the Hadamard- and the 7r/2-phase gate. The one-qubit rotations and the CNOT-gate 
alone form a universal set of gates. The reason why we also include the Hadamard- and the ir/2- 
phase gate is that here they are treated differently from the general rotations, as can be seen from 
the propagation relations (O) and (§3). By propagation through a Hadamard- or 7r/2-phase gate, the 
gate is left unchanged while the byproduct operator changes; whereas for the propagation through a 
general rotation, the rotation changes and the byproduct operator stays the same. Thus, for finding 
the byproduct images the general rotations in Af can be replaced by the identity. Only the CNOT-, 
Hadamard and 7r/2-phase gates have an effect. The special treatment of the Hadamard and the ir/2- 
phasc gate is advantageous with respect to the temporal complexity of a computation, because if one 
uses the propagation relation ( |T^ ) the implementation of the Hadamard- and the 7r/2-phase gate does 
not need to wait for results of any previous measurements. To sum up, to each possible J\f belongs a 
unitary operation U(Af) in the Clifford group and a corresponding matrix C(Af), such that 



2 (C(AT)I) = U(N) 1 (I) U(Afy , VI e V. 



(31) 



Let us now give the propagation matrices for propagation through CNOT-, Hadamard and 7r/2-phase 
gates. The propagation matrices C are conveniently written in block form 







\ 


Cxx 


Czx 




C X z 










) 



c 



\ 

where C xx , C zx , C xz and C zz are n x n matrices with binary- valued entries. 
For the Hadamard gate on the logical qubit i one finds 

[C ra (ffW)] H = [C„(ffW)] H = 6 kl + 5 ki S a , 
[C zx (H^)] kl = [C xz (H^)] kl = S ki S a , 



(32) 



(33) 



where e.g. \C XX (HW) \ u denotes the entry of row k and column I in C xx . Note that the qubit index i is 
not summed over in ( pS] ) and that the addition is modulo 2. 

For the 7r/2-phase gate £/J (tv/2) on the logical qubit i one finds 



C xx (U z i) (tt/2)) 

C' xz (uP(t:/2)) 
C zx (U^(w/2)) 



= Sm , 

= Ski, 

= SkiSu, 

= 0. 



(34) 
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For the CNOT-gate on control qubit c and target qubit t one finds the propagation matrix C(CNOT(c, t)) 
with 

[C xx (CNOT(c,t))] kl = Ski + S kt S ol , 
[C zz (CNOT(c,t))] kl = 6 kl + S kc 5 tl , 

CUCNOT(c,t)) = 0, 

C xz (CNOT(c,t)) = 0. 

We will make use of the propagation matrices in the discussion of temporal complexity of algorithms on 



the QC C in Section 5.2 



For the action of the propagation matrices C on the vectors I € V there exist conserved quantities. 
One of them, I^ 1 I Z: 2 + \Ix,2, is discussed in the next section. 

3.4.4 Conservation of the symplectic scalar product 

The symplectic scalar product 

(Il,I 2 ) S = I f 1 - ' I 7 .;T,-.2 mod 2 (36) 
is conserved. For any Ii , I 2 £ V and C the identity 

(Ii,I 2 ) s = (CTi,CI 3 ) 5 (37) 

holds. Let us briefly explain why the symplectic scalar product ([36]) is conserved. First, note that 
the symplectic scalar product tells whether two operators 1(J-i), 2(h) in the Pauli group commute or 
anti-commute, 

I(Ii)I(Ia) = (-l) (Il,l2)s 2:(I 2 )2:(Ii). (38) 



Relation ( p8[ ) is the only place in this paper where we pay attention to the sign factor of a byproduct 
operator. There, the product, e.g. I(Ii)Z(l2), denotes the usual operator product. However, everywhere 
else in this paper a product Z(Ii)Z(l2) denotes operator multiplication modulo a global phase factor ±1, 
i.e. the product is normal ordered as in (f26[ ) and the phase factor is dropped. 

Using relation (fj8[), the invariance (|37|) of the scalar product (|6|) is easily demonstrated. Consider 
the quantity X(CIi)J(CI 2 ) with I(CIi) = Ul{Ii) W and I(C*I 2 ) = U2(h) W as in (|TJ). Then, we find 

I(CIi)J(CI 2 ) = U 2(h) U^U 2(h) 

= Ul(h)I(l 2 )W 

= (-l) (Ilj2) ^C/ 2(h) 2(h) W (39) 

= (~l) (Iul2)s U2(h)U^U2(h)W 

= (-lf lJ ^2(Ch)2(Ch), 

where the third line holds by (|38|). On the other hand, as we can see from (|38| ) directly that 

2(Ch)2(Cl 2 ) = (-l) (CluCl2)s 2(Ch)2(Ch)- (40) 



From ( |39[ ) and ( |40P together it follows that (Ii,I 2 ) s = (CIi,Cl2) iS , as stated in (|37 

The symplectic scalar product (J3TJ) will prove useful in determining the measurement angles from 
previously obtained measurement results. 

3.4.5 The cone test 

The cone test is used to find out whether two measurements, which are part of some gates of a circuit, 
influence each other, i.e. whether one of the measurements has to wait for the result of the other. The 
cone test does not reveal which of the two measurements has to be performed first. 
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Let j, k be some cluster qubits k £ C and j G . Qubit j is not measured in the first measurement 
round and thus the observable measured at qubit j is a nontrivial linear combination of a x and cr y , hence 
j can be in the forward and backward cones of some other cluster qubits. We would like to find out 
whether j is in the forward or backward cone of k. For this question the cone test provides a necessary 
and sufficient criterion. It reads 

VkeCJeQ™: j g fc(fc) v j g bc(fc) (F^F^s = l. (41) 

To check whether a qubit lies in some other qubits backward or forward cone we only need the two 
forward images and can use the symplectic scalar product. 
We further observe that 

V j, k G Q (1) : j € fc(fc) k G bc(j). (42) 

If we confine fc to fc G C C wc can insert ( fl2| ) into ( f4l| ) such that the expression on the l.h.s. of (|d|) 
becomes symmetric with respect to j and k. This fits in well since the r.h.s of (041) is also symmetric. 

Similar to ( p4| ) we can give a criterion for whether or not a qubit j G is in the forward- or 

backward cone fc(<?), bc(<?) of some gate g. It reads 

V.g G JV, j G Q (1) : j G fc(.g) V j G bc(.g) (Fj,F ff )s = 1. (43) 

The proofs of the cone tests (|4l]), (fl|]) and the forward-backward cone relation ( f42| ) are given in 
appendix |A|. 

3.5 To what a quantum logic network condenses 

Simulating a quantum logic network on a QC C is a two-stage process. Before the genuine computation, 
we feed a classical computer with the network to be simulated. It returns the quantities needed to run the 
respective algorithm on the QC C . These quantities are the sets Qt of simultaneously measurable qubits, 
the measurement bases of the qubits k G Qq, the algorithm angles V 3 / algo f° r ^ G C\Qq, the backward 
cones bc(fc) of the qubits k G Qq, the byproduct images Fj for j G C and the initialization value li n it of 
the information flow vector I(t). Together these quantities represent the program for the QC C . 

In |l(| we wrote that the set of one-qubit measurements on a cluster state represents the program. 
Now we can be more specific about the measurement pattern representing the program for the QC C . The 
measurement pattern has both a temporal and a spatial structure. The temporal structure is given by 
the sets Qt of simultaneously measured qubits. The spatial structure consists of the bases (cr x -, a y - or 
a z -) of the measurements in the first round and of the measurement angles in the subsequent rounds. 
The measurement angles can be determined only run-time, since they involve the random outcomes of 
previous measurements. The measurement angles are determined using the algorithm angles and the 
byproduct images. 

4 Computational model for the QC C 

In the preceding two sections we have established the notions of the sets of simultaneously measurable 
qubits, backward cones, byproduct images, measurement angles and the information vector. In this 
section, the computational model underlying the QC C is described in these terms. First, wc would like 
to give a summary of the characteristic features of the model: 

• The QC C has no quantum input and no quantum output. 

• For any given quantum algorithm, the cluster C is divided into disjoint subsets Q t C C of qubits, 
t = 0, 1, . . . , imax, where Q p n Q q = for p ^ q and Ut=cT Qt = C. These subsets are measured one 
after the other in the order given by the index t. In measurement round t the set Qt of qubits is 
measured. 
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Figure 2: General scheme of the quantum computer via one-qubit measurements. The 
sets Qt of lattice qubits are measured one after the other. The results of earlier mea- 
surements determine the measurement bases of later ones. All classical information 
from the measurement results needed to steer the QC C is contained in the information 
flow vector I(t). After the last measurement round i ma x, I (Wax) contains the result of 
the computation. 



• The classical information gained by the measurements is processed within a flow scheme. The flow 
quantity is a classical 2n-component binary vector where n is the number of logical qubits of 
a corresponding quantum logic network and t the number of the measurement round. 

• This vector I(i), the information flow vector, is updated after every measurement round. That is, 
after the one-qubit measurements of all qubits of a set Qt have been performed simultaneously. 
I(£ — 1) is updated to through the results of these measurements. In turn, l(t) determines 
which one-qubit observables are to be measured of the qubits of the set Qt+i- 

• The result of the computation is given by the information flow vector I(t m ax) after the last mea- 
surement round. From this quantity the readout measurement result of the quantum register in the 
corresponding quantum logic network can be read off directly without further processing. 

We should make a comment on the first point. The QC C has no quantum output. Of course, the final 
result of any computation -including quantum computations- is a classical number, but for the quantum 
logic network the state of the output register before the readout measurements plays a distinguished role. 
For the QC C this is not the case, there are just cluster qubits measured in a certain order and basis. If, to 
perform a particular algorithm on the QC C , a quantum logic network is implemented on a cluster state 
there is a subset of cluster qubits which play the role of the output register. These qubits are, however, 
not the final qubits to be measured, but among the first (!). 

The QC C has no quantum input. This means that the quantum input state is known and can thus 
be created from some standard quantum state, e.g. |00...0), by a circuit preceding the main part of 
the computation. Shor's algorithm where one starts with an input state l/\/2(|0)i + is an 

example for such a situation. Other scenarios are conceivable, e.g. where an unknown quantum input is 
processed and the classical result of the computation is retransmitted to the sender of the input state; or 
the unmeasured network output register state is retransmitted. These scenarios would lead only to slight 
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modifications in the computational model. They are, however, not in the focus of this paper. The reader 
who is interested in how to read in and process an unknown quantum state with the QC C is referred to 




4.1 Obtaining the computational result from the measurement outcomes 



Now that we have defined the information vector I in (|25|) and have seen that the result of the computation 
can be directly read off from the x-part of I, we will explain how I depends on the measurement outcomes 
{sfc} and the set {/t/c} of binary numbers that determine the cluster state \4>)c m For this purpose, 
we will express Us R in terms of {sk} and {k^}, and use the isomorphism ( p6| ) to obtain I. As in (p4|), 
Usr can be decomposed into two parts, the accumulated byproduct operator Us and the contribution 
from the "readout" measurements, U R . We will proceed in three steps. First, we will consider for 
the case where only the cluster qubits for the "readout" and those of which observables cos<p a x + simp a y 
are measured are present. Second, we will extend the obtained result to the case where both the relevant 
and the redundant qubits are present, i.e. where a universal cluster C is used for computation. Third, 
we will include the contribution U R from the "readout" qubits. 

To derive I as a function of {sfc} and {ftfc}, we need to define the following sets. C is a universal 
cluster. Let O C C be the subset of the cluster which, in the simulation of a quantum logic network on 
the QC C , consists of the readout qubits. Let Cn C C denote the cluster that contains only the relevant 
cluster qubits, i.e. those which are measured in a direction in the equator of the Bloch sphere, and the 
"readout" qubits. Be Q ,z C C the set of qubits of which the operator a z is measured. Among these sets, 
the following relations hold: 

Cn U Qo,z = C 

c N nq^ z = o. 

Further, let us denote a standard form of the cluster on which the gate g can be implemented as C(g). 
The cluster C(g) shall consist only of essential qubits, i.e. those which are measured in a direction in the 
equator of the Bloch sphere (x-y-plane). The set Cn of qubits is the union of all clusters C(g) 

Cn=\J C(g). (45) 

Each C(g) is subdivided into an input zone Ci(g), an output zone Co{g) and the body of the gate Cm (<?); 
Ci{g) n Co{g) = 0,Ca/(<7) = C(g)\ (Cr(<?) UCo(g))- In the two cases described in Section |^, the general 
rotation and the CNOT-gate, these sets are C(U Rot (£,r), ()) = {1,2,3,4,5}, Ci(U Rot (^V,0) = {1}, 
C M (U Rot (C,V,0) = {2,3,4}, Co(U Rot (Z,r), 0) - {5} and C{GNOT) = {1, .., 15}, C T [CNOT) = {1,9}, 
C M {CNOT) = {2, .., 6, 8, 10, .., 14}, C (CNOT) = {7, 15}, with the labeling of qubit sites according to 
Fig.§ 

Let us now explain in greater detail the notation Kk,i ,k>i,o which made use of in equations (^), (^|) 
and ( |l(]| ) in Section ^. Some of the sets C(g) have an overlap. If a qubit k in the output zone Co (.9), 
k G Co(g) is not a readout qubit, it is also in the input zone of some gate g succeeding g, i.e. k £ Cj(g). 
The procedures to implement the gates g, g on clusters C(<?), C(g) depend on the eigenvalues of the states 
\(j>)c(g)! |</>)c((?) m the respective eigenvalue equations of form (|2|). Namely, the procedures depend on the 
set {n' a ,a G C(g)} or {n' a ,a G C(g)}, respectively. In the described case, we have k G C(g), k G C(g) 
and k G Cn- Thus, there are three eigenvalue equations associated with k, one for the state \4>)c(g), 
one for \4>)c(g) and one for \4>)c N - ^ n au the three cases the symbols denote different operators, 

because the sets ngbh(k) are different. For C(g), k G Co(g) has a left but no right neighbor. For C(g), 
k G Cj(g) has a right but no left neighbor. For Cn, k has both a left and right neighbor. We specify 
the eigenvalue in the relevant equation (^) for the state \4>)c(g)-> where k is an output qubit, by n' k Q , i.e. 

K^\4>)c{g) = (— l) f " fc - \4>)c(g) ■ Similarly, we specify the eigenvalue in the equation for \4>)c(g), where k is 
an input qubit, by n' k l and in the equation for the state \4>)c N by n' k . The n' k 7 , n' k and n' k are generally 
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different, but related in the following way 



Wk\3g,g £ TV, such that fc G Co{g) A k e Ci{g) : K k = K k I + n' k Q , 

Vfc I 3g G TV, s.th. fc G Co (g) A ^3.g G TV, s.th. k G d(g) : n' k = n' k \ (46) 

Vfc | -Bg G TV, s.th. fc G Cote) A 3 .9 e ^ s - th - fc e C K<?) : 4- = «4j- 

The latter case of (|46| ) applies to cluster qubits which belong to the cluster equivalent of the input- or 
output quantum register. The first line of (^) is proved in appendix |^, the second and third line are 
straightforward. 

For all qubits k in the gate bodies, i.e. for which 3g G TV such that k G Cm(<7) C Cjv, the n' k that 
defines \4>}c( g ) and the n' k that defines \4>)c N are equivalent, 

Vfc | 3. 9 G TV, such that k G C M (.g) : 4ec( 9 ) = 4ec„ (47) 

and are therefore denoted by the same symbol. Provided with these definitions and relations, we can 
now start to express I in terms of {s k , fc G C} and {n k , k G C}. 

Let us -in the first step- discuss the accumulated byproduct operator Us for a computation on 
the special cluster Cjy. To Us contribute all the byproduct operators Us{g) that are created in the 
implementation of the gates g. For all necessary cases, the general rotations dq), the CNOT-gatc ( [To| ) 
and the special rotations Hadamard-gate and 7r/2-phase-gatc ([l3|), the byproduct operators Us(g) can 
be written in the form 

u*(jg)=l n (u k y k+K ' k A I n (UkT k+K 'A ( n (^r^l^o?), (48) 

V'6Cj( g ) / \feec M (g) / V fceC °(f) / 

where we have attributed the contributions to Uy, which depend on n' k j, n' k or n' k to the qubit fc. 
Uo(g) is constant in the measurement outcomes {s k ,k G C/v\0} and {t^fc G C N \0} and we therefore 
attribute it to the gate g as a whole rather than to a particular cluster qubit. For all rotations we have 
Uo(g) = 1, but for the CNOT-gate -if realized as depicted in Fig. [TJ- the contribution is nontrivial as can 
be read off from ©, U {CNOT) = . 

To determine the effect of U^(g) on Uz we propagate, by use of the propagation relations (pM|), ( p"2| ) 
and (|l~4|), the byproduct operators U^(g) forward to the cut f2 which intersects the corresponding network 
TV just before the output. The forward propagated byproduct operator that results from the byproduct 
operator U?,(g) we denote by Uy(g)\a- In the same way, the forward propagated byproduct operator 
originating from Uk, the byproduct operator generated via the measurement of qubit fc, is denoted by 
Uk\n for all fc G Cn\0 . Finally, the forward propagated byproduct operator originating from Uo(g), 
the byproduct operator attributed to the gate g as a whole, is denoted by Uo(g)\n- To give an explicit 
expression, be O the vertical cut through a network TV at the output of a gate g and U(No^q) the 
unitary operation in the Clifford group which corresponds to the part of the network TV with all the 
one-qubit rotations except for the Hadamard- and 7r/2-phase gates replaced by the identity, as explained 



in Section 3.4.3, Then, the forward propagated byproduct operators are given by 

Uv(g)\ n = U(Mo^n)Ux(g)U(Mo^V 

U k \ n = U(Afo^n)U k U(Afo^n) 1 (49) 
Uo{g)\n = U{No^n)U {) {g)U{No^ 

The contribution U^(g)\n from the gate g to Us is 

Uz(g)\n=l II u k\n k+K ' k A ( II ( II U (g)\n- (50) 
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The total byproduct operator Us is the product of all forward propagated byproduct operators Us(g)\ti, 
u ^ = I\ g£ Af U v(g)\n, and thus given by 

g&M \ fceCj(g) k£C M {g) fcSCo(fl) / 

As can be seen from eq. (|5l|), for k G Cm(<7)j 9 € A/" the accumulated byproduct operator C/g depends only 
on the combinations Sk + K k . This is, in fact, true for all k G C as will be shown below. That this property 
of Us in its general form is not directly visible in eq. (|5l|) is a disadvantage of this equation. A further 
disadvantage of ( |5l| ) is that it depends on the n' k T , k € Cj{g), k[, I € Cm (g) and n' m Q , m G Co(g), but 
not on the K k , k G C which define the cluster state \<f))c- 

To find a more convenient expression for Us we first observe that by use of eq. ( |5l| ) Us can be written 
in the form 

Us=[ ]J U *\n k -const ({4, k e C N }) . (52) 
\feeCjv\o / 

This is possible because in eq. ( J5l| ) for k G Co (.9) the factors contributing to C/s do not depend on 
{sfc}. Thus, if one only picks the Sfc-dependent contributions the product index variable k runs over 
U g eAf Ci(9) u Cm (9) = Cn\0. Now, Us depends for all k G Cn\0 only on the combinations s k + n' k 
which can be seen as follows: Let us consider two cluster states \(j>)c N and \<j>)c N on the cluster Cm 
which are related via \4>)c N = °z \4>)c N > such that the respective eigenvalues are related by iz' k {\4>)c N ) = 
K' k (\<p)c N ) + 1 mod 2, and K;(|0)c N ) = K /(l ( / ) )civ) f° r all I ^ k. Suppose one would use the cluster state 
\4>)c N f° r a computation. For each choice of k G Cn\0 qubit k is measured in some direction r on the 
equator of the Bloch sphere, i.e. the operator r ■ cr® with f ■ e z = is measured, and the measurement 
result is s k - For the resulting state one finds 

l + i-iyr-yW {k) l + (-ir + V.q( fc ) - 

2 \4>)c N = o\ ' j lw c " ■ t 53 ) 

The state on the r.h.s. of eq. (|5^) is modulo the posterior equal to the state into which one had 
projected if the cluster state \4>)c N instead of \4>}c N was used for computation and the result Sk + 1 mod 2 
was obtained. Since we are only interested in the measurement results but not the resulting quantum 
state, the posterior ai is irrelevant for the computation. Using a cluster state |^>)cjvj characterized by 
K k , and obtaining the measurement result Sk is equivalent to using a cluster state \4>)c N , characterized 
by K k + 1 mod 2, and obtaining the measurement result s k + 1 mod 2. Thus, Us can depend only on 
Sk + n' k mod 2 for all k G Cm\0, i.e. 

Us= I J] U k Sk+ < J • const. (54) 



Now, the constant operator can be identified via comparison of eqs. ( |5l| ) and (|5J). For this purpose, we 
choose {«4 j = 0, k[ = 0, n' m = 0, Wg G Af, k G Ci(g),l G Cu{g),fn & Co(ff)} and insert it into eq. (|5l|). 
This choice implies via eq. ( ^ ) that {k'^ = 0, Vfc G Cat} which we insert into eq. (|54|). In this way we 
find const = Y[g£M^o(g)\n and thus 

Us = 1] C/feS ' I+K ' fc • II ^o(ff)|n. (55) 

Let us now -in the second step- include the effect of the qubits / G Qo,z\0 on Us- These qubits 
are among the redundant qubits which are measured in the first measurement round. Redundant here 
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always means redundant with respect to a given quantum logic network to be simulated. If one starts 
with a universal cluster state \<j>)c on a cluster C and projects out the qubits I G Qo,z\0 the resulting 
state on the unmeasured qubits k G Cm is again a cluster state \4>)c N - The eigenvalues that specify the 
state \4>)c N hi an equation analogous to eq. depend on the results of the tr z -measurements on the 
qubits I G Qo,z\0- For the set \n' k , k G Cn} that specifies \4>)c N one finds 



Vfc eC N : K k = n k + ^2 s j, 

je nbgh(k)DQ ,»\O 



(56) 



as can be easily derived from eq. (|2|). In eq. (^6|) the n k are those which specify the universal cluster 
state \4>)c in (0). If one now inserts eq. (p6[) into (p5|) one obtains 



II U k \ n Sk+h 
kec N \o 



n n u k \n 

keC N \0 j\ jg nbgh(k)A 

j e Q ,z\O 



J 



n u o(9)u 



(57) 



Therein and in the following it should be understood that a product of operators is set equal to the unity 
operator if the index variable runs over an empty set. The second factor in eq. (57) can now be rewritten 
in the following way 



n n w 

k€C N \0 j\j£nbgh(k)A 
j 6 QoA° 



= II II W 

j€C N \0 k\jenbgh(k)A 
k 6 Q ,z\O 

II t^lo** 
{j-,k)\ jeC N \OA (58) 

j e nbgh(k) A 

k e Q , z \O 

ii n 

feeQo, z \0 j| j e nbgh(k)A 

j e CjvVO 

In the first line of ( p8| ) the labels j and fc were interchanged and the relation j G nbgh{k) <^> fc G nbgh(j) 
was used. In the second and third line the order of the products over k and j was interchanged. 

We now define the forward propagated byproduct operators U k \n for qubits k G Qo tZ \0 = C\Cn as 



U k \a= II U M-. 

j\ j 6 nbgh(k) A 
J G Civ\0 



VA: G Qo.AO. 



(59) 



In this way, we have traced back the forward propagated byproduct operators for qubits k G Qo,z\0 to 
those for qubits j G Cn\0 which are already known. On the level of the corresponding byproduct images 
we find via the isomorphism (E6J) 



' k = 



E 



Fj , Vfc G Q ,z\°- 



(60) 



If we insert the definition 



j\ j e nbgh(k) A 

into (p8) we obtain 



n n u„\ a «= n 

A;eCjv\0 j\ j enbgh(k)A keQ 0:i \O 

j e Qo,AO 



(61) 
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Thus, with Qo,z\0 = C\Cn from cq. (Q) and substituting ( |6l| ) into (|57|), one finds 

ux=l n u k \n sk ) I n ^in Kfc ) f n^)i« • ( 62 ) 

Finally -in the third step- we investigate the contribution from Ur to U-^r which comes from the set 
O of "output" qubits. With eq. we have Ur = Y\l where the n cluster qubits in the set O 

shall be labeled in the same way as the n logical qubits. Then, we can define the propagated byproduct 
operators Uk\n for k G O as 

U k \ n =4 fc) , VfceO, (63) 
and the corresponding byproduct images via ( ^7|) 

F fe = (0i,0 2 ,..,0 fc _i,l fc ,0 fc+ i,..,0„;0„+i,..,0 2 „) T , Vfc e O, (64) 



where e.g. lfc denotes a 1 at the fcth position of F^, k G 0. Combining cqs. (62) and ( |63| ) in (|24j), the 
extended byproduct operator Usr, = UsUr becomes 

fc(EC\Q ,, / \ge^ J 




Via the isomorphism ( |26| ) and using the definitions of the byproduct images (|27| ) and ( |28|) one can now 
express the information vector I = Z~ 1 (Usr) as a function of the measurement results {s^, k € C} and 
the {Kfc, k G C} defining the cluster state \4>)c on the cluster C, 

I = Y, s kF k + ^kFk+^Fg (66) 

To derive the expression (|66| ) for the information vector has been the primary purpose of this section. 

With the expression ( p6[) at hand we are finally able to define the quantity which carries the algorithmic 
information during the computational process and which has already been mentioned on earlier occasions 
in this paper, the information flow vector I(t). 

Definition 5 The information flow vector l(t) is given by 

i(t)= skFk+ K k F k + Y F 9 ( 67 ) 

The quantity l(t) is similar to I as given in (^6|), but to I(t) only contribute the byproduct images of 
t 

qubits from a subset (J Qi of C. The information flow vector I(i m ax) after the final measurement round 
£max equals the information vector I, 

I = I(W)- (68) 

As will be shown later, during all steps of the computation, except for after the final one, the information 
flow vector determines the measurement bases for the cluster qubits that are to be measured in the next 
round. After the final round it contains the result of the computation. Thus, it has a meaning in every 
step of the computation. No further information obtained from the measurements is needed. In this 
sense, the information flow vector can be regarded as the carrier of the algorithmic information on the 

QC c fl 

1 The way we use the term "algorithmic information" has nothing to do with the - in g eneral non-computable- algorithmic 
information content of an object as it is defined in Kolmogorov complexity theory (l4|. 
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The information flow vector has a constant part which does not depend on the measurement results 
{sfe}. This part alone forms its initialization value linit, 



E 

fceC\Q ,* 



such that I(t) becomes 



E ] 



I(*)= linit + E SfcFfc - 

fce U Qi 



(69) 



(70) 



From eq. ( |69| ) we see that the byproduct images of the gates F g do not form an independent part of the 
information specifying a quantum algorithm on the QC C . Instead, they arc absorbed into the initialization 
value Ii n it of I(<). 

The measurement bases in which the results are obtained -referred to implicitly in ( p6[ ) and (|67|)- 
are not fixed a priori, but must be determined during the computa tion . T hey will be calculated using 
the byproduct images {F&, k £ C} and I(t), as explained in Sections L2 and O. Besides the byproduct 
images, the algorithm angles (/Sj.aigoi 3 £ Q' 1 ' are a l so needed to determine the appropriate measurement 
bases. They are related to the network angles ^qinj 3 S that specify the one-qubit rotations in the 
corresponding quantum logic network via 



¥>j,algo = (-1) <y3j,qln . 3 G Q 



(1) 



where r)j is given by 



The pair of equations (0), ( 
become apparent in Sections 



E 



fe| k e C\Q ,z 

j e bc(fc) 



E 

g\ a^Af, 

j S be (g) 



1. 



(71) 
(72) 



^2h is, for the moment, just a definition of the algorithm angles. It will 
.2 and 4.3 that this definition is indeed useful. 



4.2 Description of the model 



As already listed in Section a quantum algorithm on the QC C is specified by the sets Q t of simul- 
taneously measured qubits, the backward cones bc(fc) of the qubits k £ Q , the measurement bases of 
the qubits k e Q , the byproduct images Fj for j e C, the algorithm angles <^ ia i go for I G and 
the initialization value I; n it of the information flow vector l(t). If an algorithm is not given in this form 
but rather as a quantum logic network composed of CNOT-gates and one-qubit rotations, the above 
quantities can be derived from the network as explained in the previous sections. 

Let us summarize this step of classical pre-processing. First, the measurement pattern is obtained -if 
one has no better idea- by patching together the measurement patterns for the individual gates displayed 
in Fig. 0. This gives the measurement directions for the qubits k £ Qq. The network angles (fj q i n for 
the qubits j £ are taken from the quantum logic network to be simulated. To determine the sets 
{Qt, t = 0..i max }, we need the forward cones. The forward cones fc(fc) for all qubits k £ C can be obtained 
using the expressions (||), (|l(]) for the byproduct operators and the propagation relations ([ll]), (|l^) and 
(|i~4|). From the forward cones we derive a strict partial ordering "-<" (|l^) among the cluster qubits, and 
from the strict partial ordering we derive the sets Qt C C via (|2^). The byproduct images Ffe for the 
qubits k £ C\Qo,z are obtained from their definition ( p7j ) once the corresponding forward propagated 
byproduct operators are obtained from (f49| ). The byproduct images of the qubits k £ Qo,z\0 are traced 
back to those in the set C\Qo,z via e q- (p0[). The byproduct images of the remaining qubits in C, k £ O, 
are given by ([34j). To determine the algorithm angles we need the backward cones bc(fc) for the qubits 
k £ Q and the backward cones of gates bc(g). Then, the algorithm angles are given by (|7l|), (|72|). 
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Finally, for the initialization value Ij n ; t of the information flow vector we need the byproduct images F g 
of the gates g which we obtain from eq. (|2q), Ij n jt is set via d6S|). All the pre-processing required to 
extract the listed quantities from a quantum logic network can be performed efficiently on a classical 
computer, see Appendix |(]. 

The scheme of quantum computation on the QC C comprises several measurement rounds in which 
the following steps have to be performed: 

1. First measurement round. 

(a) Measure all qubits k € Qq. Obtain measurement results {sk\k G Qo}- 

(b) Modify the angles y^aigo f° r the continuous gates 

¥>j,algo ► ^j.algo = ^algo (-1)^5 ( 73 ) 

with 

rfi= E s * (74) 
fceQob'ebc(fe) 

for all j G QW. 

(c) Update the information flow vector from Ij n i t to 1(0) 

1(0) = I init + E s kFk- (75) 



2. Subsequent measurement rounds. 

Perform the following three steps ( pa| ) - (p3) for all qubit sets Q t C C\Qo in ascending order, 
beginning with Q\. In the measurement round t, 

(a) Determine the measurement bases for j G Qt according to 

^mea3 = ^ ia l g o(-l) (I(t " 1),Pi)s (76) 

(b) Perform the measurements on the qubits j G Qt. Thereby obtain the measurement results 
{ Sj G {0,l}|iGQ t }. 

(c) Update the information flow vector I 

I(*) = I(*-1)+ E *i F i- ( 77 ) 

ieOt 

The information flow vector I(t m ax) after the final measurement round t mSuX equals the information vector 
I, as can be seen from (|6^). At the end of the computation, from I we can directly read off the result l x 
of the computation. I x is identical to the readout of the corresponding quantum logic network. 

Remark 1. Note that in the first measurement round the byproduct operators created by the measure- 
ments on qubits in Q have been propagated backwards to set the angles {y/. a j go }- There is also a scheme 
in which the byproduct operators caused by the measurements in the initialization round are propagated 
forward to set the modified algorithm angles {v^ajgo}- I n that scheme, the update of the information flow 
vector I(i) and the rule to determine the measurement angles v?j,meas are the same as in the described 
scheme, given by ([76|) and (|77|). What is different is the initialization and the appearance of a step of 
post-processing. In the modified scheme, in eqs. ([72]) and (f74j) the backward cones bc(fc) are replaced by 
the respective forward cones fc(fc) and I; n jt is set to zero. The quantity which was Ii n jt in (J6FJ) is computed 
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as well but now stored as an auxiliary quantity AI until the end of the computation. After the last mea- 
surement round t mSuX , the information vector I then is obtained by the relation I = I(t max ) + AI, which 
requires the extra post-processing step and extra memory during the computation. We have chosen to 
present the scheme with backward propagation of byproduct operators in order to avoid this superfluous 
post-processing. This way, the quantity I(t) which steers the computational process directly displays the 
result of the computation after the final update to I(i m ax)- 

Remark 2. This comment concerns the choice O = O of the cut on which the byproduct images 
Ffc and F g are evaluated. In the visualization of the QC C as an implementation of a quantum logic 
network the cut SI plays a distinguished role. The byproduct operators accumulated at fl determine how 
the "readout" measurements have to be interpreted. In the computational model underlying the QC C , 
however, the former readout qubits are just qubits to be measured like any other cluster qubits. Here, 
the cut 51 is not distinguished. Due to the invariance ( |37| ) of the symplectic scalar product (p6| ) the 
byproduct images Ffc, which enter the expression © for the <^fc,meas directly and via ( |69| ) and (|77|), can 
be evaluated with respect to any vertical cut O through the corresponding quantum logic network. The 
information vector I which displays the result of the computation in its x-part 1^ would then be related to 
the information flow vector after the final measurement round I(t m ax) via I = C(No^n) I(^max)- Thus, 
the particular vertical cut O = H was chosen just to avoid an additional step of post-processing. The 
dependence on the cut O would vanish altogether if one would write the n output bits of the quantum 
computation in the form [I x ] i = (I|e>, fi\o)s f° r suitably chosen {/) € V,i = e.g. for the case 

= n, f x = (0,..,0;1,0,..,0) T , f 2 = (0, ..,0; 0,1,0, ..,0) T , and the other f h i < n accordingly. 



4.3 Proof of the model 



In this section it is shown that if we run the QC C according to the scheme described in Section 4.2, we 
obtain the same result as in the corresponding quantum logic network. This requires to prove that (a) one 
does indeed choose all the measurement angles correctly and (b) obtains at the end of the computation 
the result 1^, the x-part of the information vector I as given in (|66|). 

To show point (b), we use (|69|), ( |75| ) and (fn]) and obtain for the information vector 

i = y K ^F k + y f q + s kF k + Y skFk 

fceC\Q . 2 g^M fc6Qo fee|J'^ ax Q; 

which coincides with (|66|). This ensures that we obtain the right vector I at the end of the computation, 
provided the measurement bases were chosen appropriately, as required for (a). This is checked below. 
First we observe that the measurement angle y?j 5 meas 

and the network angle <Pj A i n are for all j e 

related in the following way 

tineas = (-1) ' ^,qln , (78) 

with 

= y, s k + H K k + Y 1 mod 2 - ( 79 ) 

kec\jefc(k) fceeXQo,* |jefc(fc) se.V|jefc(g) 



Why docs the pair of equations (T7q), (79) hold? As can be seen from the propagation relation for rotations 



Til), the network and the measurement angle of a qubit j € QW can differ only by a sign factor ±1 and 



can therefore always be related as in (78). The first and the third sum in (|79|) follow from the definition of 



the forward cones of the cluster qubits and of gates in Section 3.2. The measurement angle at j acquires 
a factor (— l) Sk if j £ fc(fc) and a factor of (—1) for each gate g with j £ fc((?). The K-dependent part 
can be derived in the same way as it was derived for J7e through equations ([32]) - (^2|). Note that $j 
depends, similar to the information vector I, only on the combination Sk + Hk for k £ C\Qo, z and on Sk 
for k £ Qo,z- 
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Now, we rewrite the quantity iDj in the following way 

■dj = ^ Sk + E Kk + E 1 mod 2 

fcec|jefc(fc) fe6C\Qo,= b'6fc(fc) geA/"|jGfc(g) 

E s * + E K fe + E 1 + 

feeC|iefc(fe) fcec\Qo.z|jefc(fc) ge.AA|jefc(g) 

E s k+ E Kfc + E 1 mod 2 

vfe6Qob'ebc(fc) feeC\Q ,= b'6bc(fe) ge-VliGbcfg) / 

E + E Sfc + E Sfc + ( 80 ) 

feeQob'efc(fc)vjebc(fe) keQW \jek(k) fceQolj'Gbc(fc) 



Si S2 S3 

+ E ^ + E Kk + 

fcec\Q . 2 b'efc(fe) vjebc(fc) fcec\Q .z b'ebc(fc) 

V * ' V v ' 

S4 S5 

+ ^2 1 + E 1 mod 2 - 

g6J^|jefc(s)vjebc( ff ) gGA^b-Gbc(s) 



5s S7 

We now discuss the seven terms Si, ■■ ,Sj. All sums are evaluated modulo 2. 
Term Si of @: 

Si= E Sfe = E s fc (F fc ,Fj)s, (81) 

fee<3ob'efc(fe)vjebc(fe) fceQc 

where the last identity holds by the cone test (|4l|). 
Term S 2 of ©: 

Let be j £ Qt and k e Qi. Qubit j can only then be in the forward cone of k, j £ fc(fc), if i < t. Hence 

s 2 = e Sk 
= E 

t-i 

fce|J Qiliefc(fc) 

= ]T Sfc (F fc ,F,)s- (82) 

t— 1 

fee U Q; 

In (|2|) the last line again follows by using the cone test 
Term S 3 of ©: 

S 3 = E s k = v'y (83) 

feeQob'Gbc(k) 



This equity follows by the definition of 77J in (|74|). Thus, the term S3 is the contribution to i9j coming from 
the first measurement round where the algorithm angles {c/^.algo} are changed to the modified algorithm 
angles {^- algo }. 
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Term S* 4 of ph 



^4= Y Kk= J2 ^k(F k ,Fj) s , (84) 

feec\Q ,= | jefc(fe) viebc(fe) fc6C\Q ,z 



which follows by the cone test ([5|). 
Terms S 5 + S 7 of (fcj): 



s 5 + s-= y Kk + S 1 = 7 &- (85) 

fceC\Q ,*|j£bc(fc) 9GAT|jebc(g) 



via the definition ( [72] ) of the 77^. 
Finally, term S 6 of @: 

S 6 = Y, 1 =Y( F 9> F As> (86) 

ffGA/"|jefc( ff )vjebc(g) seJV 

which follows by the cone criterion ( f43| ) for gates. 

Now we combine these seven terms Si, ..,87. By (|8(]) - ( |S6"|) we obtain 

^ = + + E ( F s> F ^s + Y MFfc.F.Os + £ s fe( F ^ F j)s 

9£Af fceC\Q t-i 

fce U Qi 

= % +^ + (I(i-l),F,) s . (87) 

The last line follows from the definition ( |67| ) of the information flow vector. If we consider the relations 
©' (III) and © between the angles (f jA i go , ^, a ig and P^meas, we find 

^rncas H-l)^-M algo . (88) 

Now we insert (|87]) into (^8|) and obtain 

</?j,meas — ^j,algo V L ) > 

which proofs that the assignment of the measurement angles (fq) is correct, and thereby concludes the 



proof of the computational model described in Section 4.2 



5 Logical depth and temporal complexity 

The logical depth has, to our knowledge, only been defined in the context of quantum logic networks, 
but it can straightforwardly be generalized to the QC C . In networks one groups gates which can be 
performed in parallel to layers. The logical depth of a quantum logic network then is the minimum 
number of its layers. Similarly in case of the QC C , one can group the cluster qubits which can be 
measured simultaneously to sets Qt- There, the logical depth of the QC c -realization of an algorithm is 
the minimal number of such sets. 

Since the one-qubit measurements on the cluster state mutually commute, one may be led to think 
that they can always be performed all in parallel. They could, but then the measurements would in 
general not drive a deterministic computation. 

In the following, we will denote the logical depth in the context of the QC C by D and the logical 
depth of a quantum logic network by D^j-. 
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5.1 D = 1 for circuits in the Clifford group 

The Clifford group of gates is generated by the CNOT-gates, the Hadamard-gates and the 7r/2-phase 
shifts. In this section it is proved that the logical depth of such circuits is D = 1 on the QC C , independent 
of the number of logical qubits n. For a subgroup of the Clifford group, the group generated by the CNOT- 
and Hadamard gates alone we can compare the result to the best known upper bound for quantum logic 
networks, where the logical depth D_\f scales like O(logn) |l3| ]. 

The elementary gates we use are the Hadamard gate H = \j\[2{a x + a z ), the 7r/2-phase gate 
U z (n/2) = cxp(— in/4:<T z ), and the CNOT-gate between neighbouring qubits, whose realization on the 
QC C is depicted in Fig. |l|. Out of the latter we construct the CNOT-gate between arbitrary qubits via 
the swap-gate composed of three CNOT-gates. Hence, out of the elements displayed in Fig. [l| any circuit 
in the Clifford group can be composed. At this point we must emphasize that in a practical realization of 
a QC C we would not perform a general CNOT in the described manner using the swap gate. There is a 
more efficient realization for the general CNOT, whose spatial resources scale more favourably This gate 
will be displayed elsewhere. The purpose here is to keep the argument compact rather than the gates. 

It is possible to measure all qubits at once. This works since, as shown in Fig. [l], all cluster qubits 
necessary for the realization of the CNOT-, Hadamard- and 7r/2-phase gates are measured either in the 
eigenbasis of a x or of a y . The redundant qubits are measured in a z , as explained in |Tof] ■ Thus, none of the 
qubits is measured in a basis whose proper adjustment requires classical information from measurement 
results at other qubit sites. This concludes the proof of D = 1 for circuits in the Clifford group. In a 
computation, of course all the measurement results obtained must be interpreted. Therefore, there exists 
a contribution to the computation time from classical post-processing. The connection between logical 



depth and computation time is discussed in section 5.3 



5.2 D = 2 for circuits of CNOT-gates and a £/(l)® n -subgroup of rotations 

In this section we prove that the logical depth D of a circuit composed of either CNOT-gates and rotations 
about the x-axis or of CNOT-gates and rotations about the z-axis is D = 2. This set of circuits contains 
all circuits of diagonal 2-qubit gates as a special case. For circuits of diagonal 2-qubit gates we can 
compare our result D = 2 to the best known result []l3f for quantum logic networks where the logical 
depth scales logarithmically in the number of gates. 

Here we give the proof for circuits of CNOT-gates and rotations about the z-axis U z (a) = e~ la ^~ . 
The elementary gates used are (a) the rotations about the z-axis U z (a) = e~ ta ^~ , and (b) the CNOT- 
gate between neighbouring logical qubits. The realization of the rotation U z is depicted in Fig. |l|. Of the 
CNOT-gate between neighbouring qubits we construct the swap-gate between neighbouring qubits and 



by that the general CNOT-gate, as in section 5.1. The strategy to implement the circuit is then: (1) to 
measure all those qubits on C which are to be measured in the eigenbases of a x , a z or <x y ; and (2) to 
measure the remaining qubits, i.e. the ones which are measured in a direction in the x — y-plane. 

The result that the measurements in step (1) can be performed in one step has already been shown in 



section 5.1, It remains to be shown that the measurements in the tilted measurement directions of step 
(2) can also be performed in parallel. Let j and I be two cluster qubits which arc measured in a tilted 
basis in step 2 in order to implement the rotations. Using ( p0| ) one finds 

D > 2 3 j, I G Q (1) : I -< j (that is, Q< 2 ) ^ 0). (89) 

Further holds 

I -< j => 3 k G Q (1) : j G fc(fc), (90) 

because the strict partial ordering is generated by the forward cones, i.e. I -< j <^=> cither j G 
fc(0, or 3 (fti, ...k r ): ki G fc(i) A {k s G fc(fe s _i)| 2 < s < r} A j G fc(fc r ). 
Moreover, from the criterion (^) one derives 

jGfc(*)=>(F il F fc )s = l. (91) 
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Now, by putting the implications fl89|), ( p0| ) and (|9l|) together we obtain 

D>2^3j,keQ^ : (F i ,F fc )s = l, 



(92) 



which we negate to obtain 

Vj, k G Q (1) : (Fj, F fc )s = D < 2. (93) 

Next it is proved that (F_,-, Ffc)s = docs indeed hold for all j, k £ 

A measurement of a qubit at site k , which is part of the implementation of a rotation about the 
z-axis (central qubit 3 in Fig. p]c), generates a byproduct operator {Uk) Sk = (<7 z ) Sk . This can be seen 
from equations ([|), (^) and (g)] Note that in Fig. [l]c. qubits 1,2,4 arc measured in the cr x -cigenbasis, 
they belong to the set Qq. Now let be i the number of the logical qubit on which the rotation U z (ipk) is 
performed by the measurement of cluster qubit k. Further, let O be a vertical cut through the network 
simulated by the QC C . O intersects each qubit line only once. In particular, it shall intersect the qubit 
line i just at the output side of the rotation U x (<pk)- Thus, the image Fk\o of Uk on the cut O is 



•k\o 



( 



\ 



\ 



'kz\0 



with Fi 



k zA 



(94) 



/ 



What we see from ( |94| ) is that Fk x \o = 0. Be No^n the part of the network M which is located between 
the two cuts O and f2. The byproduct image Fj. corresponding to Uk is then given by 



fen 



C(Afo-.n)F k \ . 



(95) 



The only gates that contribute to C(Afo->n) are the CNOT gates, as described in section 3.4.3 . The 
propagation matrices for CNOT gates (|3^) have block-diagonal form. Hence, using ([50]) the propagation 
matrix for the network Afo^n has block-diagonal form 
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(96) 



From (|94[) , (|95| ) and (|96|) it follows that the x-part of the byproduct image vector Ffc vanishes for all k 



x \k 



= Vfc G Q 



(i) 



(97) 



Hence by the definition of the symplectic scalar product (|36|), we obtain (Fj,Fk)s = for all j, k G 
This proves via J93| ) D < 2. The measurements to implement the one-qubit rotations can thus all be 
performed at the same time. In (^) the case D = 1 can be easily be excluded for all interesting cases such 
that only D = 2 remains. This concludes the proof of D = 2 for circuits of CNOT-gates and rotations of 
the form e~ 41/3 ~^. The proof for circuits of CNOT-gates and rotations e~ lv ^~ runs analogously. Now let 
us discuss the special case of circuits composed of diagonal two-qubit gates. A diagonal gate Gd in the 
computational basis is of the form 



Gd 



\ 



1 / 



(98) 
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Figure 3: Network for a diagonal gate composed of rotations U z and CNOT-gates. 



modulo a possible global phase which is not relevant. 

The network of rotations about the z-axis and of a CNOT-gate shown in Fig. || realizes a general 
diagonal two-qubit gate. In order to obtain the angles ipi, ip2 and (ps specifying the diagonal gate Gd in 
(|9§|), one chooses the following angles for the three z-rotations in this network 

1 / \ 
a = -^{-^1-^2 + ^3), 

P = \ (-^1 + ^2-^3), (99) 

7 = 2 (~ lpi +ip ' 2 + (f>3 >- 

Thus, a circuit of diagonal two-qubit gates can also be regarded as a circuit of z-rotations and CNOT- 
gates. Therefore we find D = 2 for circuits of diagonal two-qubit gates on the QC C . This result can be 
compared to the best known upper bound Jl3| for quantum logic networks where the logical depth is of 
the order O(\ogno) with uq the number of two-qubit gates. 

5.3 The logical depth D is a good measure for temporal complexity 

In this section, we will express the computation time as a function of the logical depth. The computational 



model described in section 4.2 consists of an alternating series of measurement rounds and classical 
processing of the thereby obtained measurement results. The classical processing contributes to the 
duration of the computation and will therefore enter into the relation between the computation time and 
the logical depth. For the computation time, this results in a correction logarithmic in the number n of 
logical qubits involved, and thus the computation time is no longer the logical depth times a constant. 
For all practical purposes, however, this logarithmic correction is small compared to the time required 
for the genuine quantum part of the computation, consisting of the measurements. 

Let Aq be the time required to perform the simultaneous measurements in one measurement round 
and Ad the time required for the elementary steps of classical processing: say, addition modulo 2 or 
multiplication of two bits. The time T c i (t) required for classical processing after each measurement round 
has two contributions. First, the time T c ij(t) to update the information flow vector and second, the 
time T c i ±{t) to determine the signs of the measurement angles of all measurements in the next round. 
The total computation time T comp is given by 

D-l 

^comp = 

DA Q + TciM + T ci At) (ioo) 



The update of the information vector l(t) according to (|77|) can be done for all 2n components in parallel. 
The update I(t — 1) — > I(t) following measurement round t requires the time that it takes to add up 
HQtll bits modulo 2. As the drawing below illustrates, T c n(t) is logarithmic in ||Qt||. 
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The number of qubits in the set Qt is bounded from above by ||C|j since Qt C C. Here, C is any cluster 
sufficiently large to carry the network to be simulated. Thus 

T ciiI (i) < A cl log ||C||. (101) 

To determine the proper measurement angle (fik,mcas f° r the measurement on qubit k £ Qt+i in the 
next measurement round requires, according to ((76j); the evaluation of the symplcctic scalar product 
(I(t),F k )s- This requires 1 step for multiplication and log2n steps for addition modulo 2. Thus, 

T clt± = A c/ (logn + 2). (102) 

Combining ( |l00| ), ( |l0l| ) and (|102| ), the total computation time T comp is bounded from above by 

Tcomp < D A Q (l + ^ [log ||C|| + log 7i + 2]^ . (103) 

We see that, although the computation time T com p is linear in the logical depth D, it contains con- 
tributions logarithmic in the number n of logical qubits and in the cluster size ||C||. These logarithmic 
contributions are, however, suppressed by the ratio between the characteristic time for classical processing 
and the characteristic time for the von-Neumann measurements, A c //Aq. This ratio can, in practice, be 
very small. Therefore, the logarithmic corrections become important only in the limit of large clusters and 
large n. As will be argued below, even in the regimes where a quantum computer is believed to become 
useful, say n ss 10 5 , the logarithmic corrections have only a minor influence on the total computation 
time. 

We now eliminate the dependence of the total computation time on the cluster size ||C||. For this we 
assume that on the QC C we simulate a quantum logic network with the network logical depth D^/. Now, 
we give an upper bound on ||C|| as a function of n and D^f. As displayed in Fig. [l], a single CNOT gate 
has height 3 and width 6 on the cluster C. Here we do not count the output qubits of the gates since they 
also form the input qubits of the gates in the next slice. As in Fig. [|, the rotation has height 1 and width 
4, if the output qubit is again not counted for the width. The wires for the logical qubits on the cluster 
can be arranged with distance 2. Each set of parallelized gates will at most require a slice of width 6 on 
the cluster. The circuit as a whole requires an additional slice of width 1 for the output. A swap-gate 
that is composed of three CNOT-gates, requires an array of 3 x 18 qubits on the cluster. If a general 
CNOT-gate on the cluster were composed of a next-neighbour CNOT-gate and swap-gates (in practice 
it is not), then it would require at most an array of [height] x [width] =[2n — 1] x [18(n — 2) + 6] qubits. 
Hence, a CNOT gate would, to leading order, consume at most 36n 2 cluster qubits. Each rotation would 
require at most -in the worst case where on the network it could not be performed in parallel with other 
gates- a slice of width 4 on the cluster, so it consumes, to leading order, at most 8n cluster qubits. The 
total number of gates in the network is bounded from above by nD^f. The simulation of each gate costs 
at most max(36n 2 ,8n) = 36n 2 cluster qubits. Hence, the size of the required cluster is bounded by 

||C|| < 36n 3 Av. (104) 

If we now use the assumption about a good quantum algorithm that the logical depth scales polynomially 
in the number of qubits n, 

D H = cn p , (105) 
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and insert (104) and ( 



1051 ) into we obtain 



1+^(4 + log 9c) 



^(p + 4)logn 



From a practical point of view, we find that the logarithmic corrections -even for numbers n of logical 
qubits in the range of 10 5 - play a minor role since they are suppressed by the ratio A c //Aq. Wc 
could plug in some typical numbers, say Aq = 1 /is, A c ; = 1 ns, p = 3 and n = 10 5 , to obtain 
A cJ /A Q (p + 4) log ti w 0.12 (or w 0.24 for n = 10 10 ). 

The spatial overhead ||C|| is polynomial in the number n of logical qubits. But, if one adopts this 
more practical viewpoint one may not be satisfied by the mere result that the spatial overhead scales 
polynomially, but might want to know what the scaling power actually is. Above we found that ||C|| 
scales with the (p + 3)th power of n. However, in the above argument, we focused on the computation 
time where the precise value of exponent for the spatial scaling did not play an important role, thus 
have been extremely wasteful with spatial resources. A more careful discussion yields a more favourable 
scaling of the spatial overhead. 



From a strict scaling point of view, we find in (106) that the computation time is no longer equal to 



the logical depth D times a constant, but there are logn -corrections due to the classical processing. This 
is, as the above numbers illustrate, of little relevance for practical purposes. The classical processing can 
be parallelized to such a degree that it increases the total computation time only marginally. 



6 Discussion 

The discussion about the logical depth of certain algorithms with the QC C in Section || showed that 
there are ways of information processing possible with the QC C which cannot be explained adequately 
in network model terms. This made a computational model appropriately describing the QC C desirable. 
The computational model underlying the QC C that we found does not seem to have much in common 
with the network model. It is based on objects of a different sort which require an interpretation. In this 
section, we attempt to clarify the role of the binary valued information flow vector l(t) and that of the 
stepwise measured quantum state. 

What is the role of the information flow vector I(i)? In every computational step except the final one 
the information flow vector I(i) is completely random. So one might ask whether it contains information 
at all. It does, since in every step except the last one it tells what has to be done next. After the final 
computational step at time £ ma x the quantity I(t m ax) contains the result of the computation. Thus, the 
quantity I(i) has a meaning in every computational step. In this sense, it represents the algorithmic 
information in the described scheme of quantum computation. 

What is the role of the stepwise measured quantum state? To see that explicitly, let us consider 
the scenario where a quantum computation is halted in the middle and continued at a later time by 
another person who only knows which steps of the computation are left to perform but does not know 
what has been done so far. In analogy to a teleportation protocol where both the result from the Bell 
measurement and the quantum state at the receivers side are required to reconstruct the initial state, 
the halted computation can be successfully completed only if both pieces -the intermediate information 
flow vector I(U) and the half-measured quantum state- are stored until the computation proceeds. Thus, 
the quantum state cannot be neglected just because it does not appear in the formal description of 
the computational model. The quantum correlations in the stepwise measured state are what basically 
enables the described way of quantum information processing. However, the role of this state is a passive 
one. It serves as a resource that is used up during the course of computation. 

For the description of the QC C there is no necessity to refer to the "qubit" as a basic notion of 
quantum information theory. As described above, the stepwise modified quantum state attains the role 
of a consumable resource but not that of the carrier of algorithmic information. It is thus sufficiently 
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described in terms of standard quantum mechanics, namely as a quantum state of entangled two- level 
systems. The processed algorithmic information is classical. 

Let us at the end of this discussion come back to the role which the randomness of the individual 
measurement results plays for the QC C . It may surprise that a set of classical binary numbers represents 
the algorithmic information in a scheme of quantum computation. In the network model the quantum 
state (of the quantum register) is usually considered to represent the processed information. For the 
QC C , the situation is different. There, the randomness of the individual measurement results makes it 
necessary to store classical steering information. The need to process this information has called for a 
novel information carrying quantity. What, in a network-like description of the QC C , has been regarded 
as a mere byproduct turns out to be the central quantity of information processing with the QC C . 

7 Conclusion 

We have described the computational model underlying the one-way quantum computer, which is very 
different from the quantum logic network model. 

The formal description of the QC C is based on primitive quantities of which the most important are 
the sets Q t C C of cluster qubits defining the temporal ordering of measurements on the cluster state, 
and the binary valued information flow vector I(t) which is the carrier of the algorithmic information. 

The information processing with the QC C goes beyond mere emulation of quantum gates by sequences 
of measurements. The complete description of the computational process on the QC C contains the tem- 
poral order in which the measurements are performed, and the most efficient temporal order does not 
follow from the rules that apply for the temporal ordering of gates in the network model. Thus, the 
network picture is insufficient to describe computation with the QC C and is therefore abandoned. In the 
proposed scheme the unitary gates from some universal set can not be taken for the elementary building 
blocks of quantum computation. 

As a practical implication of the new model, the logical depth of certain algorithms on the QC C is 
lower than for their corresponding network realizations. As shown for circuits of diagonal gates or of 
CNOT-, Hadamard- and 7r/2-phase gates, the logical depths for the QC C are constant in the number of 
gates and logical qubits, namely one and two. The best bounds for networks that have been known so far 
scale logarithmically. It therefore seems that, at least for the QC C , the question of temporal complexity 
must be revisited. The tools for this discussion together with first results have been provided in this 
paper. 
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A Proof: cone test and forward- backward cone relation 



Here we prove the cone test ( |4l| ) and the forward-backward cone relation ([42|). Considering the cone 
test, first note that whether a one-qubit rotation at some position in the network is about the z-axis or 
a about the a;-axis can be identified by the potential byproduct operator produced when the rotation 
is implemented. This can be seen by inspecting ([?]), (R), (^J) and the procedure to implement a general 
rotation as described in section ||. The a;-rotations U x (£) and U X (Q of ?/r(£,77, C) m (3) are implemented 
by measurements on the qubits 2 and 4 of a 5-qubit chain. As can be seen from (8), they contribute 
to the byproduct operator Us of the rotation Ur with a x S2+Si where S2 and S4 are the results of the 
measurements on qubits 2 and 4. Further, the rotation about the z-axis, U z (£), is implemented by 
measurement of qubit 3. The contribution to the byproduct operator which is thereby generated is, from 
(@), <J z Sa . We see that x-rotations only generate byproduct operators o x and z-rotations only generate 
byproduct operators cr z . 

A byproduct operator generated via the measurement on the cluster qubit k must be propagated 
either forward or backward to possibly reach the rotation on the logical qubit i implemented via the 
measurement on the cluster qubit j. Let be Ok and Oj two cuts through the network which intersect 
each logical qubit line only once. More specifically, Ok intersects the qubit line i just before the rotation 
implemented by the measurement at cluster qubit k. Oj intersects the qubit line i just before the rotation 
implemented by the measurement at cluster qubit j. 

There are two cases which can occur. Either the cut Ok is before the cut Oj in the network Af which 
we denote by Ok < Oj, or Oj is before the cut Ok which we denote by Oj < Ok- It can also be that 
both is true at the same time but it cannot be that neither of the two relations hold. 

Case l:0 K <0.j. 

The byproduct operator generated via the measurement at qubit k must be propagated forward to 
possibly affect the measurement at qubit j. It is not possible that the result of the measurement on qubit 
j has an effect on the measurement basis chosen at k. 

Let us introduce a further cut Oji which is the same as Oj, except for that it intersects the line of the 
logical qubit i in the network Af just after the rotation implemented via the measurement on the cluster 
qubit j. The potential byproduct operator which is generated via the measurement on cluster qubit k 
and then propagated forward to the cuts Oj and Oj>, is denoted by Uk\o., and Uk\o jn respectively (the 
byproduct operators which are actually generated are {Uk\o.,) Sk and (C4|e> 7 ,) Sfc ). Further, we denote the 
restriction of the byproduct operators Uk\o.j and Uk\oj, to the logical qubit i by [f/fc|oj]j and [f/fe|c> J ,] .. 
The two cuts differ only on the logical qubit i, and there only by the side of the respective cut on which 
the rotation is located. Therefore, using (|Tl"|), it follows that Uk\o. r — Uk\o ■ Hence also 



Pi 



k Oj 



[U k \c 



(107) 



If the rotation implemented via the measurement on cluster qubit j is about the a;-axis, then the measure- 
ment on qubit j has to wait for the measurement on cluster qubit k iff [J/fc|c>jL contains a contribution 
cr z . The measurement on j itself produces a potential byproduct operator [?7j|©j/]. = & x - Similarly, if 
the rotation implemented via the measurement on j is about the z-axis then the measurement on j has 
to wait for the measurement on k iff [C/fc|o J ] i contains a contribution <r x . The measurement on j itself 
produces a potential byproduct operator [C^/|oj,] ■ = cr z . 

Because of (107) Uk can as well be evaluated at the cut Oj> instead of Oj. The byproduct operator 
on the intersection of qubit line i and cut Oj> resulting from the measurement on qubit j can be written 
in the form 



j\Oj> 



with 




for z-rotations 
for x-rotations 



(108) 
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The byproduct operator on the intersection of qubit line i and cut Oji resulting from the measurement 
on qubit k reads 

[^l^] i =(^ ) ) XM (^ ) )* k,< - (1 ° 9) 
One can now easily check for both the cases of an x- and a z-rotation implemented by the measurement 
on qubit j that the measurement of qubit j must wait for the result of the measurement of qubit k iff 

Xj,iZk,i + Zj,iXk,i = 1 mod 2. (110) 

Now we note that the correspondence between ( X -'' 1 ~\ and [fjlo,/]^ an d between f Xk ' 1 ~\ and 

[ffeloj,] . is via the restriction of the isomorphism (|26| ) on qubit i. Thus, Xjj, are just the i-components 
of Ixlow and l z \Oji, respectively. Equivalent relations hold for Xk,i,Zk,i- One finds 

Xj,i=[lx,j\o J ,] i , Zj,i=[lz,j\o JI ] i _ / ni x 
%k,i = [ix.klOj^ , Zk,i — [i«,Jfe|Oj/] j 

Further we observe that 

[IxjM^O, [I,, i |o J ,] i =0forain^i, (112) 

since the byproduct operator introduced by the implementation of the rotation acts, at the cut Oji, 
non-trivially only on the logical qubit i. Thus we can write 



•Kj.iZk.i ~t~ Zj^Xj~ i — ^ ^ XjjZfc i -\- ZjjXfcJ 

= (t\ ,Mo,) s ( 113 ) 
= (Ij|n,Ifc|n) s 
= (Fj,F fc ) s , 

where the second line holds by the de finit ion ( |36|) an d the third by (^) and the conservation (^) of the 
symplectic scalar product. Inserting (113) into ( |110| ) yields 

O k <Oj: jefc(t)^(Fj,F k ) s = l. (114) 

For Ok < Oj, j G bc(fe) cannot occur, hence with ( |ll4|) , 

Ok < Oj : j G fc(fc) V j G bc(fc) «=> (F J; F fc ) s = 1. (115) 

Case II: Oj < K - 

First we observe that j can only be in the backward cone of fc, but not in the forward cone. Thus, the 
byproduct operator generated via the measurement on k must be propagated backwards in the network 
to reach the gate for whose implementation qubit j is to be measured. The reasoning is completely 
analogous to case I, up to the fact that the potential byproduct operator generated via the measurement 
of cluster qubit k is in this case propagated backwards onto the cut Oj> . Qubit j is in the backward cone 
of qubit k iff the quantity (lj\ojn^-k\Oj,) s is equal to 1. Again, by conservation ([37]) of the symplectic 
scalar product follows 

Oj<0 K - jebc(fc)^(F,-,F fe ) s = l. (116) 



For Oj < Ok, j € fc(fc) cannot occur, and therefore with (116), 

Oj < O k : j G fc(fc) V j G bc(fc) <=> (F is F fc ) s = 1. (117) 



Now we combine the two cases and with (115) and (117) we obtain 

k G C, j G Q (1) : j G fc(fc) V j G bc(fc) <=> (F,-,F fe ) s 
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which proves the cone test (|4l|). 

The proof of the cone test for gates ([l3]) goes along the same lines, only the byproduct operator (U k ) Sk 
generated via the measurement at cluster qubit k G C has to be replaced with the byproduct operator 
Uo(g) of the gate g. 

Finally, the proof the forward-backward cone relation shall be outlined. Suppose that j € fc(fc). With 
the same methods as used in the proof of ( [H| ) one can derive that 

j G fc(fc) (I k \ j,Ij\oj)s = 1, /, is s 

fcebcO') (l fc | OK) l J -K) s = i. 1 8j 



Then, with ( |118[ ) and the invariance (|37]) of the symplectic scalar product 

j G fc(fc) <J=^ fc G bc(j), 

which proves (^2|). 



B Relation among K r kI , tz' k and k\ 



k 

Let us briefly explain why the first line of eq. ( p| ) holds. Be the left neighbor of the qubit k G Co(<?)nC/(g) 
in question denoted by A; — 1 and the right neighbor by k + 1. The interaction S;v takes a product state 
I-P)cjv = ®ieCiv I"!")' *° a c l us t er state |^)cjv- Since all pair-interactions in the Hamiltonian generating 
Sn commute, Sn can be written in the form Sn = Sfe-i.fc Sk.k+i S', where Sk-i.k describes the pair- 
interaction between k— 1 and k, Sk,k+i the pair-interaction between k and fe + 1, and S' all the remaining 
pair-interactions. One can easily check that 



q S k ) of _ I -WPAk-l) (k) 

q „( k ) at _ ( 



<3fc-l,fc Oz o, 



7 (fc) (fe+1) 

(119) 

(fe+1) „Y _ (fe+1) 



fe-l,fe 



with a,j9,7 G {0, 1}. Since the state \P)c N obeys the eigenvalue equation \P)c N = &x \P)c N the state 
\4>)c N obeys \4>)c N — Sn cr^ S N \4>)c N - Thus, via cqs. (||) and (11S) we can identify a = n' k . Further 
holds 

S N Sj, = S fc _i,fe5fe,fe + i Sj+^-i = (-l^?- 1 ^*? +1 > , (120) 
such that a = n' k = (3 + j. For the cluster state \<t>)c{g) on the cluster C(g) where only the qubits k—1 and 



k are present, but not qubit k+ 1, there holds |0)c( ff ) = S k -i yk a^ S\_ l k \4>)c( g )- Via cqs. @j and (119) 



we can thus identify (3 = K ko - In the same way, if we look at the cluster state \<p)c(g') on the cluster 
C(g') we can identify 7 — n' k 7 . Combining the relations for a, f3, 7, we finally obtain n' k — n' k j + n' k as 
in eq. (ph. 



C Temporal complexity of the classical pre-processing 

The time that it takes for a classical computer (i.e. a compiler) to translate an algorithm into a machine- 
specific set of operations (i.e. the machine code) is usually not regarded as to count for the temporal 
complexity of that algorithm. For quantum logic networks this viewpoint is certainly justified because 
there the complexity to compute the circuit layout is well understood and known not to exceed the 
complexity of the quantum logic network itself. 

For the QC C however, the situation is less clear. For the time being, we do not know of any other 
way to obtain the quantities characterizing a QC c -computation than to derive them from the network 
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formulation of the respective algorithm. Therefore, we must exclude the possibility that for the QC C the 
algorithmic complexity of a quantum computation is shuffled from the genuine quantum part of the 
computation to the classical pre-processing, and that this classical pre-processing may be exponentially 
hard. As will be shown below, such a case does not occur. All the classical pre-processing can be done 
in polynomial time. 

To see this, we assume that the quantum algorithm on n logical qubits is given as a sequence of 
||7V|| elementary gates. For good quantum algorithms, \\Af\\ is polynomial in n, as is ||C||, the number of 
physical qubits in the cluster C required to run the algorithm (see Section |5.3| ). 

The layout of the measurement pattern requires to assign ||Qo|| measurement bases and ||C\Qo| 
angles. Creating the pattern is for itself not a problem since it can be obtained by patching together 
the measurement patterns of the elementary gates which are available in block form. The temporal 
complexity for this step is thus 0(j|Cj|). 

To obtain the byproduct images we introduce \\Af\\ vertical cuts Oi, i = 1, .., \\Af\\ to the network, one 
after each gate (such that 0\\at\\ = ^) an d compute the 2n x 2n-matrices C(.A/© 4 ->n) for i — 1, .., \\Af\\ — 1, 
starting with i = \\Af\\ — 1. The operational effort for this is of the order O (tt. 3 1 1 _A/"| | ) . By use of these 
matrices the byproduct images for cluster qubits k £ C\Qo, z can now be obtained via (p^), which requires 
0(n 2 ) elementary operations per byproduct image. The way to obtain the byproduct images F g of the 
gates is the same. For k £ Qo. z \0 at most four byproduct images have to be added in (|60|), which 
requires 0(n) operations. The computation of for k £ O is trivial. Thus, to compute a byproduct 
image requires at most 0(n 2 ) operations per cluster qubit or gate such that the complexity to compute 
all of them is at most O (n 2 (\\C\\ + \\N\\)). 

The backward- and forward cones of the cluster qubits k £ C are computed using the temporal ordering 
of gates in a sequence representing the quantum logic network and the cone test (ffl|). The number of 
cone tests that have to be performed in each case is \\Af\\ (||A/"|| — l)/2 where the computational effort for 
each test scales like 0{n). Thus, the complexity of this step is 0(n||A/"|| 2 ). 

The forward cones generate the anti-reflexive semi ordering "-<". The semi ordering can be computed 
from them in 0(J|Cj| 5 ) steps. 

For each set Qt there have to be HQ*-*- 1 !! < ||C|| test of the relation j -< k, j £ performed to check 
whether some qubit k £ is in Q t . Also, HQ^H qubits have to be checked for each Q t . At most ||C|| 
sets Q t exist such that the operational effort to obtain the these sets is 0(||C|| 3 ). 

As far as the stated upper bounds are conclusive, it looks as if the computation of the anti-reflexive 
semi ordering is the toughest part. However, as elementary a relation between the cluster qubits "-<" is, 
for the conversion of a quantum logic network into a QC c -algorithm it needs not be computed. Please 



note that the semi ordering is finally only needed to compute the sets Qt via (20). But instead of 
computing "-<" from the forward cones and the sets Qt from "-<", the sets Qt can also be computed from 
the forward cones directly. For this, please note that 3j £ Qt\ j -< k £ Q^ <^=> 3j' £ Qt\ k £ ic(j'). The 
direction "<^ = " is obvious with j = j' . The opposite direction, "=>", holds by an argument analogous 
to the one justifying (|9^). In fact, statement "=>'' is the same as <M) with QW replaced by QW. Thus, 
cq. ( p0| ) can be replaced by 



Qt = {k £ Q (t) \ -a? £ Q (t) : k £ fc(/)} . (121) 



To set the algorithm angles via (|7l|), ( |72[) requires at most ||C|| + ||7V|| additions per angle and there are 
at most ||C|| such angles. Hence, in total it takes O (||C||(||C|| + ||A/"||)) operations to set them. Finally, to 
initialize the information flow vector via ( |69|) requires 0(n(||C|| + ||A/"||)) operations. Thus we see that 
all classical processing requires only a polynomial overhead of elementary operations and can therefore 
be done in polynomial time. 
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